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Abstract 

This paper gives an exposition of well known results on vector parti- 
tion functions. The exposition is based on works of M. Brion, A. Szenes 
and M. Vergne and is geared toward explicit computer realizations. In 
particular, the paper presents two algorithms for computing the vector 
partition function with respect to a finite set of vectors / as a quasipoly- 
nomial over a finite set of pointed polyhedral cones. We use the developed 
techniques to relate a result of P. Tumarkin and A. Felikson (and present 
an independent proof in the particular case of finite-dimensional root sys- 
tems) to give bounds for the periods of the Kostant partition functions of 
Eq, E7, Eg, F4, G2 (the periods are divisors of respectively 6, 12, 60, 12, 
6). 

The first of the described algorithms has been realized and is publicly 
available under the Library General Public License v3.0 at 
|http: //vectorpartition. sourcef orge .net/ i We include (non- unique) 
partial fraction decompositions for the generating functions of the Kostant 
partition function for A2, A3, A4, B2, B3, C2, C3, G2 in the appendix. 

1 Introduction and notation 

Given a finite set / of non-zero integral vectors with nonnegative coor- 
dinates (with respect to some basis), and a vector 7 with coordinates 
(7^, . . . , 7"), the vector partition function Piij) is by definition the num- 
ber of ways we can split 7 as an integral sum with non-negative coeffi- 
cients of the vectors in /. It is well known that there exist finitely many 
pointed (i.e. with walls passing through the origin) polyhedral cones, 
with walls parallel to hyperplanes spanned by subsets of I of rank n — 1, 
such that Pi is a quasipolynomial over each as a function of the coor- 
dinates (7^,..., 7"). We will reprove these results (CoroUarv I4.3|l : for 
original proofs, history, and further information on the subject we direct 
the reader to (C PV OS] , |Stu95) and the references therein. 

Our definition of quasipolynomial (Section [3TT} is done via quasinum- 
bers (=linear combinations of indicator functions of affine translates of 
lattices) and is slightly different from the commonly accepted one. The 
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commonly accepted definition of quasipolynomials uses holomorphic func- 
tions of the form e^'^'^"*'' and we relate it to ours in the remarks after 
Definition 13 . II and in the appendix (indicator functions are not holomor- 
phic) . 

In Section |3] we present an elementary algorithm for computing the 
vector partition function as a quasipolynomial. The algorithm allows us 
to prove that the partition function is a quasipolynomial over a finite set 
of closed pointed polyhedral cones (Theorem 14. 2p . but doesn't prove that 
the walls of the closed pointed polyhedral cones are parallel to hyperplanes 
spanned by subset of rank n — 1 of the starting set I. 

In Sections[2]and[3]we present another algorithm, based on the Szenes- 
Vergne formula ( SV031 Lemma 1.8]). The algorithm is implied in |SV03| 
but is not formulated in the exact same fashion as here; the current expo- 
sition is geared towards explicit computer realization. A similar approach 
is taken in [CP05| . The walls of the chambers of quasipolinomiality pro- 
duced by this algorithm are parallel to the hyperplanes spanned by subsets 
of / of rank n— 1, but a priori may fail to pass through the origin. However, 
two quasipolynomials coincide if they take the same value over an affine 
translate of a pointed polyhedral cone of non-zero n-volume (Lemma l3.ip 
and so the combinatorial chambers produced by the algorithm in Sections 
[2]and[3]are "glued" by those produced by the Section |4] algorithm and the 
other way round. This proves Corollary 113] (see [BBCVOep . In practice, 
the algorithm in Sections [2] and O has turned out much faster than the 
one in Section |4l 

The algorithm in Sections [2] and [3] is based on the notion of a partial 
fraction decomposition in the ring Q[a;i, . . . , Xn, . . . , -^] [ jr^lagz" 
(Definition 12. 3p . The terminology "partial fraction decomposition" ap- 
pears in the same context in [SV03) and |CP05] . however our definition is 
different to avoid computations with complex roots of unity (which eases 
computer realization of the algorithm). In Section [2.31 we show how to 
use the elongated Szenes-Vergne formula (Lemma I2.1f c)) to produce a 
partial fraction decomposition. In Section [3] we show how to obtain the 
power series expansion of the decomposed partial fractions by computing 
explicitly the form specified by |BV99I Theorem 1]. 

The algorithm in Sections[2]and|3]allows us to give upper bounds on the 
periods of the Kostant partition functions of iSg, E7, Eg, F4,, G2 by proving 
that they are divisors of respectively 6, 12, 60, 12 and 6 (Proposition 
IS.ip . The periods for the classical root systems were first computed in 
[BBC V06] ■ and are 1 for type A and 2 for the other series. Proposition 
15. H is based on the key Lemma [5.21 which was first proved using different 
techniques in [FT07| . 

Section fS. 1.1 [ presents a modification of Section [23[ applicable to partial 
fraction decompositions of the generating function of the Kostant partition 
function of the classical root systems of type A, B, C and D (see Section 
ll.l.ip . The modified algorithm Section |6]TT] has turned out in practice to 
work considerably faster than the straightforward application of Section 

121 

We have realized the algorithms presented in Sections [2l [3] and [6] un- 
der the name "vector partition function" program. The program and its 
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source code are publicly available at |http : //vectorpartition . sourcef orge .net/l 
under the Library General Public License v3.0 (free use, modification, and 
redistribution as long as the same rights are granted to recipients in all 
consequent redistributions). Section [6.21 gives a very brief description of 
some of the features the program and the appendix presents printouts of 
(non-unique) partial fraction decompositions (as formatted by the pro- 
gram) of the generating function of the Kostant partition function for A2, 
A3, A4, B2, B-i, C2, C3, and G2. 

The current version of the "vector partition function" program (as 
of September 2009) can compute algebraic expressions and combinatorial 
chamber subdivisions for the Kostant partition function of A^, A4,, A3, A2, 
B4, B3, B2, C4, C3, C2, Di and G2. The software package LattE can also 
compute algebraic expressions for the vector partition function. LattE was 
released in 2003 by a team directed by J. A. De Loera at University of Cali- 
fornia, Davis, and is available at http: //www, math .ucdavis . edu/~ latte/ software .htm) 
Another project that computes vector partition functions is S. Verdoolaege's 
program "Barvinok" , http://www.kotnet.org/~skimo/barvinok/ For 
a detailed description of the algorithm of the "Barvinok" program see 
[BCVWOSpl . Software for computing values of the Kostant partition func- 
tion for the classical root systems has also been written by M. Vergne and 
C. Cochet in MAPLE. 

Acknowledgements. The author would like to thank Michele Vergne, 
Thomas Bliem and Alexey Petukhov for the valuable advice and discus- 
sions on vector partition functions, as well as Pavel Tumarkin for the 
valuable discussion on root systems and Coxeter groups. 

1.1 Notation 

The notation is mostly standard with the exception of the vector space/dual 
space notation. An ambient vector space of dimension n will be fixed and 
denoted as C". Greek letters will denote vectors, except for the number tt 
and the letter r which will always stand for quasinumber. For a preferred 
choice of basis ei, . . .e„, we will denote by e* the corresponding elements 
of the dual vector space defined by e*(£j) ~ 5ij- Then the expression 
refers to a;^^'''' We will use the notation 7' := e*(7) for the 

i"' coordinate with respect to the preferred basis Ei. Whenever we use 
the preceding three notations, a preferred choice of basis will be either 
explicitly stated or implied. 

Q[a:i, . . . ,Xn,^,. . . , ^][Y:r^]c«ez" will denote the subring of Q(a;i, 
. . . , Xn) generated over Q by a;^, a;"'^ and j^-^ for all a G Z". 

1.1.1 Root systems 

The Kostant partition function is the vector partition function with re- 
spect to the vectors with non-negative coordinates of sets of integral vec- 
tors called root systems, (coordinates are chosen with respect to a special 

^Up to the author's knowledge, a detailed explanation of the exact relation between the 
partial fractions approach and the Barvinok algorithm approach for computing vector partition 
functions is missing from the literature. 
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basis called "simple"). Root systems play a central role in the representa- 
tion theory of Lie algebras and classify the simple finite-dimensional Lie 
algebras. There are 4 classical series of root systems, A„, Bn, Cn, Dn, (in- 
dexed by a parameter n) and five exceptional root systems, Ea, E7, Es, F4,,G2- 
In the literature positive root systems of the root systems An, B„, Cn and 
Dn are often realized in non-simple basis coordinates as 



An,n>2 


: A+ = 


{vi - m 


j < j G {1,... 


,n + l}}; 




Bn,n>2 


: A+ = 




i < j G {1,... 


,n}}U{r?, 


i G {l,...n}}; 


C„,n > 2 


: A+ = 


{■m ± m 


j < j G {1,... 


,n}}\{0}; 




D„,n > 4 


: A+ = 




i < j G {1,... 


,n}}. 




and in simple 


basis coordinates 


as 






An ' 


£1 := rji 


- '72, • • ■ , 


£n . — ^7^ ^7^ 






Bn : 


£1 := Vi 


- ??2, . . . , 


En-l Tjn-l 


'^Jl ; £ n ■ — 


77"; 


Cn '• 


£1 ■— m 


- ??2, . . . , 


£n-l Tjn-l 


'^Jl ; £ n ■ — 


2?7n; 


Dn : 


£1 ■■= Vi 


~ '72, . . . , 


£n-l rjn-l 


5 £ n ■ — 


7)„_i 4- 



List of the positive root systems of exceptional Lie algebras Ea, Er, Eg,, 
Fa and G2 can be found at 

|http: //www . liegroups . org/dissemination/spherical/explorer/rootSystem. cgi| 

(the Atlas of Lie Groups and Representations team) and are automatically 
generated by the "vector partition function" program. 



1.1.2 List of specific notation 

^ 1 1 1 r 1 



• TM,c,d(7) 



1 if M7 - c = ( mod d) 



see Section fl.ll 

see Definition l3.1l TM,c,d 



otherwise 

is called quasinumber, the Af is a m x n integer matrix, c is an m x 1 
integer column vector. 

Bk{x) = J2t^o > 1, -^0(3;) = a; -I- 1 - Bernoulh sum polynomial 
in X. We consider Bk as an element of C[x]. 

B]^^''^{x) := '^^^QTi,m,d{t)'t'' - T-BernouUi sum, m,l £ Z>o (consid- 
ered as 1 by 1 matrices), d G Z>o, a; G Z. 

&ff^''^{x) := X]t=o '''i,m,d{'t)'t'^ - the extension of the r-BernouUi sum 
over a:: G Q via the right hand side of l|16p . where m, I G Z>o, d G Z>o. 

. See Section Fl.ll 



• a:: := Xi . . . : 

• a — 

axj 

• Aij({Qi}fli) = {aiQfi -I- • ■ • -|- ajvQivlai G -R} - cone (respectively 
lattice) spanned by the a^'s over R, where R is an arbitrary subset 
(resp. a subring) of C. 

• Pai,...,an • C" — >■ Z - vector partition function with respect to the 
vectors Oi. 
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• Z" := Az({ei}) - short notation for the lattice generated by the 
vectors Ei, where Si is the imphed preferred choice of basis. 

• (•,•) - bihnear form defined in Section[3l 

• - symmetric bihnear positive definite form on a root system 
induced from the Kilhng form and normahzed as in the proof of 
Lemma l5.2l 



2 Partial fraction decomposition algorithm 
2.1 Partial fraction decomposition definition 

Definition 2.1 A closed pointed polyhedral cone C is a subset of Q" 
of non-zero n-volume given by homogeneous non-strict rational linear in- 
equalities. The hyperplanes given by the defining inequalities are called 
supporting hyperplanes of C on the condition that they intersect C m a 
set of non-zero n— 1-volume. The intersection of a supporting hyperplane 
with C will be called a wall of C. 

Remark. A wall of a closed pointed polyhedral cone is in turn a closed 
pointed polyhedral cone in its corresponding supporting hyperplane. 

Definition 2.2 A fraction indexing set is defined to be a finite collection 
of triples of s -tuples {(ai, . . . , Qs, li, . . . ,ls, mi, . . . , rris)}, where s £ Z>o, 
all cii £ Z", all li G Z>o, all rUi G Z>o. 

The following definition is similar to |CP05I Proposition 3.2]. 

Definition 2.3 Let q £ Q[xi, ... ,x„, ^, ^][j^]asz" ■ A fraction 
indexing set for which all participating s-tuples ai, . . . , of vectors are 
linearly independent is called a partial fraction indexing set. A partial 
fraction decomposition of q indexed by a partial fraction indexing set I 
is defined to be a finite set of polynomials Vi = {Pai,...,as,ii,-.-i3,m.i,...,ms 
\ {ai, . . . ,as,li,. . . ,ls,mi, . . . ,ms) £ 1} C Q[xi, . . . ,Xn, ^ , . . . , ^] such 
that 

ll,...,ls, rni,...,ms) 
mi,...,ms)£J 

The set of pointed polyhedral cones {Aq^^ ({ai , . . . , Qs}) | (ai,...,Qs, 
li, . . . ,ls, mi, . . . , TTis) e / for some m^, li} will be called cone support of 
the partial fraction decomposition Vi. The finite set of all vectors par- 
ticipating as generators of the elements of the cone support of Vi will be 
called support of Pj. 

Remarli. We note that the "elongations" U in the above definition could 
be substituted for a single "greatest elongation" by the geometric series 
sum formula ("the elongation formula"). However, this would surely yield 
great computer memory loss in realizations of the algorithm in Section[2]3] 
below. 
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2.2 Formulas used 

Lemma 2.1 

(a) (Szenes-Vergne formula) Let ai,...ak £ such that "^^cxi =^ 
Then 

1 1 



1 - 1 - X"!' 



(i_a:<^i)...(i_a:<^j-i) (1 - a;"3+i ) . . . (1 - ) 
("ftj ("elongation formula" (geometric series sum formula)) 

1 - a;" ~ 1 - X"" 

/or n e Z>o, and 

1 -a;-" a"" 



1 - a;" 1 - a;"" 

/or n € Z<o. 

('c^ ('t/ie elongated Szenes-Vergne formula) Let k > 1 and ai,...ak £ 
Z", a; € Z^^o sMc/i t/iat OiOi ^ T/ien 

1 1 _ 

1 - a;«i " ' 1 - a;"fc ~ 

1 /t^ ...a;"^-i"^-i PaAx'^') 



l_a;Eai«i (l-a;"i)...(l-a;"^-i) {1 - x°'i+^) . . . {1 - x"'') 

where pn(x") are the numerators in the elongation formula (h). Re- 
mark. When k = 1 (c) degenerates to (b). 

(d) 

111 



1 - a;" 1 - a;P 1 - x"-p \1-xp 1-x 



1 x"-P+x' 



_l_ 



1 - a;« 1 - a;'3 1 - x"-'^P \ l-x^ 1 - ; 

(e) Letp:= -x°'{x-^ + ■■■ +x-'^'^). Then 
1 1 



(1 - a;")' {l-xPy 



' 1 
V - 



3^ (1 — a;")* (1 — a;"-"'5)'+'"-* V m — 1 



Z + m - t - 1 



^ (1 - a;/3)* (1 - a;«-"'')'+'"-* V / - 1 
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Proof, (a) is a partial case of [SV03I Lemma 1.8]. (c) is a combination 
of (a) and (b). Formulas (d) and (e) are used in Section [S] and have been 
implemented together with formula (c) in the "vector partition function" 
program. The coefficients in (e) come from counting the number of ways 
one can get from point (t, 1) to point {l,m) in a square grid by moving 
up and right with unit steps. □ 

Remark. One could realize a formula similar to (e) for more than 
two variables. 

2.3 An algorithm for computing partial fraction 
decompositions 

We are given a finite set of integral non-zero vectors with non-negative 
coordinates A = {ai, . . . ^om} C Z" and a finite set A of multiples of 
vectors in A. Let an element q G Q[a;i, . . . , a;„, . . . , ^] [ j-r^]a,gz" be 
given in the form 

where paij,...,<:<,^,mi,...,mfc(K) G Q[xi, . . . ,a;„, ^, . . . , ^]. The algorithm 
presented here gives us a partial fraction decomposition Vi of q with 
support A (see Definition 12. 3p . The denominators of the so produced 
partial fraction decomposition will involve in general multiples of vectors 
in A lying outside of A. 

The element q will be stored in the computer's memory as a list of 
elements called fractions, corresponding to each summand in For a 
fixed fraction the number mi will be called the multiplicity of the vec- 
tor Of,;, , and {oij , . . . , oii^ } will be called the support of the fraction. The 
fractions' supports are partially ordered by inclusion. Besides the data im- 
plied by ((U, for each fraction we store a boolean flag indicating whether 
our fraction is reduced and a boolean flag indicating whether we have a 
"preferred linear dependence" . The "reduced" flag is set only if the vec- 
tors in the fraction's support are linearly independent. If the "preferred 
linear dependence" flag is set, we store an additional data structure repre- 
senting a non-zero Z-linear dependence a-^ai-^ + ■ ■ ■ + akO-i^. = aoQio with 
a distinguished index io for which ao > 0. The vector we define to be 
the gaining multiplicity vector. 

Step Initialize the starting data by setting all fractions to "non-reduced" 
and not having a "preferred linear dependence". 

Step 1 Select an arbitrary fraction from ([2]) that has a "preferred linear de- 
pendence" . In case no fraction with a "preferred linear dependence" 
exists, select an arbitrary fraction from Q that is not labeled as 
"reduced". Choose, if possible, in an arbitrary fashion a linear de- 
pendence aoOio = aiQij H [-atoti^, ao > 0, ai - ^ 0,Vj, and store 

it as a "preferred linear dependence". 
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If no choice of "preferred linear dependence" is possible (i.e. the 
support of the fraction consists of linearly independent vectors) label 
the fraction as "reduced" and go to the beginning of Step 1. 
If all fractions from are labeled as "reduced" terminate the pro- 
gram. 

Step 2 Apply the "elongated Szenes-Vergne formula" (Lemma 12. l|) ('c'l and 
the "elongation formula" (Lemma [5TlJb)) as indicated below 



apply Lcmm j2. W c) here 
/ \ "10 

1 



(l-a;"'o) 



( the rest of the multiples). 



\ apply Lcmm j2. 1 ^ bl with elongation here/ 



This procedure will erase the selected fraction from the computer's 
memory and substitute it with new list of fractions. All newly pro- 
duced fractions that have the same support as the originally selected 
fraction are set to inherit the "preferred linear dependence". The 
newly produced fractions with a strictly smaller support are labeled 
as "not having a preferred linear dependence" . 

Remark In the "vector partition function" program the application 
of Lemma [2. 11 (b) is made after all applications of the Szenes-Vergne 
formula are completed. This increases computational speed signifi- 
cantly, but requires additional accounting. 

Step 3 Add the vector aoaig to the set A. Go to Step 1 with the so modified 
list of fractions q and set A. 

Lemma 2.2 The above algorithm wtll come to a halt by labeling all frac- 
tions as reduced. The so produced sum will represent a partial fraction 
decomposition for the starting element q, whose support is a subset of A. 

Proof. If a fraction participating in ((2)) has a "preferred linear depen- 
dence" , Step 5 increases the multiplicity of the "gaining multiplicity vec- 
tor" while decreasing the multiplicity of one of the other vectors in the 
support of the fraction. Thus, after finitely many steps, the support of all 
newly produced fractions will strictly decrease. Therefore all fractions in 
the final expression will end up being reduced. □ 

Remark (by an idea of Thomas Bliem). In case there is no 
"preferred linear dependence" chosen, there is freedom in choosing the 
index ajf] in Step 1 and one can choose any strategy. One such choice 
leads to partial fraction decompositions that involve non-broken circuit 
bases in the sense of [CP04] . Suppose one wanted to compute a partial 
fraction decomposition of ^ j"^^^ , • • • , W;^ . Suppose one knew by in- 
duction that a subset S :— en, , . . . , a.i. of the vectors T := Qi, , . . . , ai„ , 



and the order of the vectors in the preferred linear dependence 
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ii < ■ ■ ■ < In, was linearly independent. Then one would add vector 
'^^k+i ^ ^^'^ check whether the span of 5* U {ctif,^i} is maximal; if 
not, one would label vector cii^^-^ in Step 1 of our algorithm as "gaining 
multiplicity" and proceed to Step 5. If the set S U {di^^^ } were of full 
rank, then one would substitute S by it and continue with the next index. 
This "simple strategy" of choosing indices in fact yields a partial fraction 
decomposition with respect to non- broken circuit bases ( [CP04] ). By def- 
inition, an n-tuple of vectors q^^ , . . . , Ofi,, , ii < • • ■ < i„ is a non-broken 
circuit base if , . . . , Qi^ and q™ are not linearly dependent for any in- 
dex m with ik < m < ik+i- We leave that to the reader to prove that 
our "simple strategy" of choosing indices in Step 1 will produce a partial 
fraction decomposition involving denominators whose corresponding vec- 
tors form non-broken circuit bases. It is important to note however, that 
the results in Section [S] are based on a different strategy on choosing the 
"gaining multiplicity" vector (namely, the vector with a smallest possible 
in absolute value coefficient in a non-trivial integral linear relation). This 
strategy has been realized in the latest version of the "vector partition 
function" program. 

For the generating function of the Kostant partition function of the 
classical root systems of type A, B, C or D (see Section ll.l.ip one can 
considerably simplify the algorithm, see Section [6. 1.11 



3 Brion-Vergne decomposition for par- 
tial fraction decompositions 

For any vector 7 define the differential operator :— '^^''Xkdk in the 
ring Q[xi, . . . , x„, ^ , . . . , ^][j^]aez'^- Let ai, . . . ,Qn be linearly in- 
dependent vectors where n is the dimension of the ambient vector space. 
Define the scalar product {•,•) by {ai,aj) = J^^a^a^. Let Pi be defined 
as the vector whose scalar product with equals Sik- We have that 

^''-(-"^)-{ r otLTw-se -Th- 



and so we obtain a Brion-Vergne decomposition, [BV99| : 



(n((%-- 



mi — 1 ///VI — a;"! 



1 1 



(3) 



(1 - a:"!)""! ■ ■ ' (1 - a;^")™" 

Let A := Az(ai, . . . , a„). Now let us explore what happens when we apply 
5,3 to a power series ^(a;) := J2-y£A P{l)x^'- 



\7eA / 7eA 
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Now let us explore what happens when we apply multiplication by X' 



/3. 



7eA 7SA 

Let J2-y£A ''"(7)2;^ be the power series expansion of j^j^_^a„f given 

by uncovering the brackets in the formal product (1 + x"^ + x^"^ + 
...)•■• (1 + a;"" + a;2°=" +...)• Then 

^ / 1 if 7 lies in Az>o({ai}) 

1 otherwise " ^ ' 



Summarizing the above formulas, we get 

n 



j:UY[(i^^ + l)ri,-S)x^ 



n mi —1 



^(7) = n n (^^^+i)-(7-^) 



n( ^^--^-^f +r-Mr(T-5). (5) 



3.1 Realizing ([5]) with quasipolynomials 

Formula (O clearly raises the question of how to realize the function r(7 — 
8). For the linear inequalities implied by "7— 5 lies in Az>o(q:i, . . . , ctn)", 
we can deal as in ,BBCV06__. First, we fix an "indicator" vector in general 
position. Then we take into consideration only the fractions whose support 
generates cones that contain our "indicator" vector. In order for this 
approach to work we would need to start with an "indicator" vector "far 
away we are from the walls", since there is a —5 shift in the argument 
of r. This naturally raises the question, how does the quasipolynomial 
expression for the vector partition function change when its argument 
is close to a wall of a cone from the cone support of a partial fraction 
decomposition. 

The answer to the this question is well known (see 'CPV08', 'Stu95]. 
[SV03] ) - the quasipolynomial expression for the vector partition function 
does not depend on the distance from the walls. We give our own proof 
(Theorem 14. 2[l in Section [J] Of course, this immediately implies that the 
algebraic expressions for the vector partition function, which are different 
for different pointed polyhedral cones, take the same values over a common 
wall of neighboring pointed polyhedral cones. 

In order to represent in the computer's memory the function r, we store 
the linear inequalities in a separate structure and introduce the helping 
function r which "forgets" the linear inequalities: 

. f 1 if 7 lies in Az(ai, . . . , a„) 

\ Q otherwise 
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/ 6ii 



blr. 



Let B be the matrix B : — 

\ bnl • • • bnn 

Then f (7 — 5) is non-zero if and only if 



\ Oil 



B 



\ 7" - 5" 



(6) 



Set Ci := "^binS" and define di to be the smallest positive integer for 
which diCi £ Z,dibik G Z, Vfc (in other words, di is the least common 
multiple of the denominators involved in the j"'-line of ([6])). Define d to 
be the least common multiple of the di-s {d will later be called "period"). 
Remarks, d is the least common multiple of the orders of the elements of 
the group Z"/Az({ai}). Thus d is a (possibly proper) divisor of the index 
of the subgroup Az{{ai}) in the group Z" (which in turn equals det(a^)). 
The period d does not change when we replace ai, Uj by ±ai, iaizLaj for 
any two indices i and j and any choice of signs (since the lattice doesn't 
change) . 

f (7 — 5) is non-zero if and only if 



dB 



( -y' 



\ 7" 



= d 



V 



( mod d). 



(7) 



where rjv/., 



The notation 







Let M — dB be integer matrix of the above system, c be the inte- 
ger vector column of the system. Define TA(c«i,...,a„),4(7^, • • • j 7") — 

TM,c,d(7S---,7") := 
1 if 7i are integers 
and the system ((7]) 
has integer solution 
otherwise 

.,a„),5 is more suggestive mathematically, but for the computer 
realization it was more convenient to use the rAf,c,(;-notation, and that is 
why we choose the latter for the rest of this the paper. We extend the 
definition of TM.c,d to M being a m x n integer matrix and c a m x 1 vec- 
tor column. This makes sense because by definition the function TM,c.d is 
supported on the a priori defined lattice Z". 

Definition 3.1 The subring of the ring of functions f : C" — > Q gener- 
ated by TM,c,d with coefficients in Q we will call the ring of quasinumbers. 
The elements TM,c,d will be called basic quasinumbers. 

The subring of the functions f : C" — >■ Q generated by the quasinumbers 
and £* we will call the ring of quasipolynomials. 

We note that relating expressions written with complex exponents and 
quasinumbers requires extra caution. One precise way of doing it is by 
saying that to every quasinumber q over the lattice Z" and for any lattice 
refinement ^Z" one can assign non-uniquely a holomorphic function / 



of the form / = 



E 



P{aa(ei 



where the 
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P{oi}'s are polynomials with rational coefficients, such that the functions 
/ and q have the same restrictions on the lattice ^Z" C C". One can do 
that by combining the one variable formula in the appendix with integral 
linear substitutions. The non-uniqueness is clear: the zero function and 
the function e^'^'^» — 1 both correspond to the zero quasinumber for the 
lattice Z", but not for the lattice ^Z". 

Making Gaussian elimination transformations mod d to the system 
([7| does not change TM,c,d- Such transformations have been realized in 
our "vector partition function" program and are used to determine when 
TM,c,d is zero and to multiply elements TM,c,d by TM',c',d' to yield an 
expression of the same type tm" .c" ,d" ■ However, the program does not 
realize all additive relations in the ring of quasi-numbers. We finish this 
subsection with a definition and an important lemma. 

Definition 3.2 A least common period for a finite set {TMi,c^.di} we call 
the least common multiple of the di-s. Let the quasipolynomial q be written 
via the generators e* and the elements of the finite set {TMi,ci,di} ■ Then 
the least common period of the set {TMi,Ci,di} is called a possible common 
period for q. The minimum over all possible common periods of q is called 
the least common period for q or simply the period of q. 

Lemma 3.1 We are using the notation from Definition \2. 1\ If two quasipoly- 
nomials qi and q2 coincide on an affine translate C' of a closed pointed 
polyhedral cone C of non-zero measure, they coincide everywhere. 

Proof. Let A be the intersection of the lattices naturally determined 
by the basic quasinumbers participating in the expressions for q\ and 
q2 ■ Then both q\ and q2 are polynomials when restricted over each lattice 
translate (5 + A as 5 varies over the representatives of elements of the group 
Z"/A. Now the statement of the lemma follows from the corresponding 
statement for polynomials. □ 

3.2 An algorithm for computing the vector par- 
tition function with respect to a set / of non-zero 
integral vectors with nonnegative entries 

• Step 1. Apply the algorithm from Section [2.31 to get a partial frac- 
tion decomposition of the generating function of the vector partition 
function, n^e/ r^- 

• Step 2. Fix a combinatorial chamber of interest and apply formula 
([5} to the partial fraction decomposition obtained in Step 1. One 
needs to fix an vector in general position which indicates the combi- 
natorial chamber. 
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4 Partition function algebraic expressions 
do not depend on the distance from the 
walls 

This section presents an independent, elementary, but computationally 
inefficient algorithm for calculating the vector partition function, whose 
aim is to prove that the vector partition function does not depend on the 
distance from the walls of the combinatorial chambers. 

Suppose we want to compute the vector partition function Pai,...,a^ 
with respect to the vectors ai, . . . , Om-i, Om and have already computed 
the partition function with respect to the first m — 1 vectors (call it 
Pai,...,a„-L-i)- Our starting point is the following completely elementary 
observation: if one chooses exactly t times the vector a„, one can de- 
compose a vector 7 in exactly P^i ,...,q„_i (7 — tQm)(7) ways. Therefore 

oc 

-Pc«i,...,Q„(7) = ^ Pai a„_i(7 - iOm). (8) 

t=0 

We will gradually refine ((51 to a complete algorithm how to compute 
the vector partition function. This elementary idea is the essence of this 
sectior0. 

Suppose as an induction hypothesis there exist finitely many closed 
pointed polyhedral cones Co, . . . , Ca/ such that Pai,...,Q„_i , restricted to 
each Ci, is a quasipolynomial; extend this quasipolynomial naturally to 
the entire space and call the corresponding element • Suppose now we 
want to compute Pai,...,a„ at a point 7 in general position, i.e. such that 
7 lies in no support hyperplane of any of the d. Let _ff be a wall of an 
arbitrary Cj such that the ray J? := {7 + tam.\t G [0, 00)} intersects H. 
By the requirement that 7 be in general position, we get that R intersects 
H transversally, i.e. if vh is a normal vector for H then {am,i^H) 7^ 0. 
Assume in addition uh is chosen to point towards the interior of Cj (i.e. 

{i'H,s} > o,yseCj). 

Definition 4.1 If {PH,am,) < (respectively {vH,otm,) > 0) we call 
such a wall H an exit [respectively entering) wall for the point 7 with 
respect to the direction am, We will use the shortened term "exit (re- 
spectively entering) wall" whenever the rest of the data is implied by the 
context. 

vh must have rational coordinates, so assume in addition that it is rescaled 
to have integer coordinates. Then the point through which 7? exits Cj is 
7 — Q; in coordinates this is the point (7"^ — oi^{^'i^v\j), . . . ,7" — 

a'"(^7V|f)). Note that the coordinates of this exit point are a linear 
rational expression of the coordinates 7'. Let Ho,Hi,...,Hn be exit 
walls and Hi, . . . ,Hn-i, be entering walls for the point 7 with respect 
to direction a„ and Co, ■ ■ ■ , Cn , N < M be pointed polyhedral cones 
such that Hi, N > i > 1 is sl shared wall of d-i and d and Hn is 

^The author has not fully understood yet the interplay between the current algorithm and 
the preceding one. 

■'The author hopes that the resulting rather cumbersome formulas will not obscure it. 
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a wall of Cn- Assume in addition that the exit wall Hm is a subset 
of a wall of Aq(q:i, . . . , cviv). Let the ray R intersect Hi in point Ai. 

Then Ai — ^ — -^—^^Ofm- Assume finally that Co U ■ • ■ U Cjv are all 

pointed polyhedral cones that have a common point with the segment 

{7 + ton.|te [0,7^ (lr^;^> Q'^]}- Set 

ti~P^^^^,l<l<N-l- to-0. (9) 

The contribution of the first chamber Co to the sum ^ is 

[tiJ L«iv-iJ Lt]v-iJ 

^Pco(7-tom)= Pcoh-tarn)- Pcoil-ta^)- (10) 

t=0 t=0 t=[tij+l 

The contribution any other chamber Cj, A*' > j > is 

[tj + lJ L*JV-lJ [tjv-ij 

Y Pc,(7-to„)= Pc,{l-ta^)- Y Pc,{l~tam). 

t=itj\ + i t=LtjJ + i t=Ltj+iJ + i 

(11) 

Note that if [tj+ij < [tjj + 1 the latter expression is zero. Summing 
over all chambers d we get the first refinement of ((8|: 

P.i......„(7) = Pc,{l)+Y\ E ^c,(7-to™)) 

3 = \( 




tam) . (12) 



Note that if the lattice spanned by the q^'s is unimodular, (for definition 
see [Stu95] l. (I12p is immediately computable via Bernoulli polynomials. 
To compute (|12|) in closed form one needs to know in advance which are 
all possible exit /entering walls for the point 7 with respect to the direction 
Qfm which in turn will account for new subdivisions of the chambers d. 

Lemma 4.1 Suppose Co is a closed pointed polyhedral eone with walls 
with rational normal vectors such that there is a unique choice of exit 
walls Ho, . . . , Hm for a point 7 in Co with respect to the direction am ■ Let 
Co, ... , Cat and Pcq, ■ ■ ■ , Pc^ '^^ above. Suppose in addition each Pc^ 
is a quasipolynomial. Then Pai,...,a,ri quasipolynomial over Co- 

Proof, (i) Assume 7 is in general position with respect to the supporting 
walls of Co, . . . , Cjv. Due to the uniqueness of the choice of exit walls such 
points 7 lie in the interior of Co- 

To prove the lemma holds for (i) all we need to do is show how to 
express (|12|) with quasinumbers. Assume the ti's are written as linear 
expressions in the 7*'s with rational coefficients represented by irreducible 
fractions and let d be the least common multiple of the denominators of 
the ti's. Assume that d > 1. 
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Let / : — >■ C be an arbitrary function. Then 



/(W) = E^i>'.<^(^)/(^-^)' (13) 

where 1 stands for the one by one identity matrix. The above is a simple 
interpolation formula: multiplying by T\^i_d(x) we guarantee that "we get 
a contribution" only if a; = /( mod d). 

We can carry out explicitly the substitution Pcii,...,a„_i|Ci (t ~ in 
the right hand side of \Vl\ using 

a-l 

TM,a.a{l - fo) = ^ T^,l ,a{t)r M ,c-^lct,a{.i) , (14) 
i = 

where 1 stands for the one by one identity matrix. Regrouping with 
respect to t (reminder: the coordinates of 7 are free variables, t is not) we 
get for each summand of (|12|l 

Ct=Uiv-iJ \ 
Ti,ra,d{t')^ | (somc quasipolynouiial in 7'). (15) 
t^LtjJ + i / 

The r-BernouUi sum B^y^''^{x) := X]t=o ^ ^ com- 
puted in closed form using quasinumbers - the d different cases based on 
the remainder of x mod d we can express algebraically with quasinum- 
bers as in (I13p . We include the formula for completeness. 



^^,m,dj-^-j ^ ^ ^ ^ ri,m,d(rd -I- s) (rd-fs) 



t = t=rdJra 
d-1 LfJ 

= Ti,^,d(s) y~^(T-d -(- s)*' 

s = r = 

a: p 

d— 1 d-1 d d 

= Ti,m,d(s) Ti.p,d(a) (rd-fs)*' 

s = p=0 r = 

d-1 d-1 

= 'ri.m,d(s) 'ri.p,d(a) 

3 = p = 

9=0 V ^ / 



where a; £ Z>o, and Bqix) are the Bernoulli polynomials with rational 
coefflcients defined by Bq(x) — '^t^g't'' G Z>o- It is to be noted that 
the numbers Ti^m,d{s) in p6|) are integer constants (zero or one) that can 
be directly substituted. Now extend the left hand side of (|16p via the the 
expression on the right hand side over x £ Q. Call the resulting extension 
^^,m,d 'pjjgjj can use our interpolation formula H13p and the properties 
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of quasinumbers to compute 



t=Ltjv-iJ 
t=LtjJ+i 

d-l 

= Tl,3,d(i]V-l)-Bj.''"''' (tjv-i) — 

3 = 

d-l 

J2^i,sAt.)B'r'' ih) ■ (17) 

3 = 

To carry out the substitutions Tj^s^ditN-i), Ti^s^ditj), we use the exphcit 
definition ^ and p4|) which completes the description of how to compute 
P2[) in closed form. 

(ii) Suppose now 7 lies on a supporting hyperplane of d for some i. 
Note that (|10[) and Hll|) continue to hold true, since if a point tj = [tj] 
belongs to more than one chamber, it will be counted only for the first with 
respect to direction Om- Therefore the algebraic expressions, obtained by 
the above procedure, will be valid for 7 lying in a supporting hyperplane 
for a chamber d, which completes the proof. □ 

Lemma [4. II leads to an elementary plan for computing the vector par- 
tition function which we summarize below. 

Step We start with a list of one polyhedral cone determined by the first 
ai, . . . , Q„ vectors (where n is rank ambient vector space). We set 
the counting index m to be n + 1. 

Step 1 We find the first polyhedral cone d over which all we have not 

computed Pa^ such that we have computed Pai....,am over 

any chamber C' that is adjacent to an entering wall for any point 
7 € Ci with respect to am- If we can't find such a chamber we 
increment the counter m. If m is greater than the largest index 
of the ai's we terminate the program, else we go to the beginning of 
Step 1. 

Step 2 We subdivide d into pointed polyhedral sub-cones C(, . . . , C; such 
that for each Cj there there is a unique choice of exit walls as required 
by Lemma [4.1l We replace the chamber d with the newly produced 
pointed polyhedral sub-cones. 

Step 3 For each newly obtained chamber C'j we compute Pqi — Pai,...,am [c'. 

in closed form via p2|) and (|4.ip . 

Step 4 We go back to Step 1. 

The fact that we can always carry out Step 2 and that the resulting 
number of chambers is finite is geometrically obvious from the fact that 
any combinatorial chamber has only finitely many faces; the detailed proof 
is left to the reader. This proves the following theorem. 

Theorem 4.2 Let ai,...,aM C Z" be non-zero integral vectors with 
nonnegative coordinates. The above algorithm will produce a finite list of 
closed pointed polyhedral cones Ci, . . . , Cm of full measure with IJilli C'j = 
Aq>o {cei, ...,q:]v) and compute the vector partition function Pai,...,a,n^Q . 
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as a quasipolynomial over each Cj. The quasipolynomial expressions are 
valid over the closures of the combinatorial chambers Cj . In particular, 
the expressions PQi,...,ajv |c '^^'^ -Pai. --,ajv |c '^'''^ equal when restricted to 
the intersection of d and Cj . 

In fact the "vector partition function" program realizes this subdivision in 
arbitrary, but small enough, dimension with arbitrary, but small enough, 
input dat£0 . An explanation of the computational difficulties of solving 
this problem for higher dimensions and higher number of faces can be 
found in |BSLV03| . 

Corollary 4.3 Let a\, . . . ,aM C Z" be non-zero integral vectors with 
nonnegative coordinates of maximal rank. 

(a) For a fixed pointed polyhedral cone d coming from the algorithm 
in the current section, the algorithm presented in Section \ 3.'A will 
produce algebraic expressions valid over the entire chamber d . 

(b) \BBCV06I Let Di, . . . , Dm be the cones spanned by n — 1 element 
subsets of a\, . . . ,aM of rank n — 1. Let d be the closed polyhedral 
cones in which the boundaries of D\, . . . Dm split Aq>o (qi, . . . Qjv). 
Then Pai,...,ajv quasipolynomial over each d. 

(c) (b) holds if we replace the cones Di by the cones in the cone support 
fDefinition \2.3]) of any partial fraction decomposition ofYl'^^i i^^a ■ 

Proof, (c) Suppose the contrary and let pj,j G J be polynomials giving 
a partial fraction decomposition for the generating function of the vec- 
tor partition function (see Definition 12. 3p given by Section 12.31 Let C 
be a fixed pointed polyhedral cone defined by the walls participating in 
the partial fraction decomposition. Let Ci^, . . . d^ be all pointed polyhe- 
dral cones given by the algorithm which have a full-measure intersection 
with C. Theorem 14.21 gives us then quasipolynomials measuring the vec- 
tor partition function over the entire d^ . On the other hand, we know 
from Section [3] that the quasipolynomials given by the algorithm in Sec- 
tion [T^l measures the vector partition function over the intersection of all 
affine translates 5 + C, where S varies over the exponents of all monomi- 
als participating in the polynomials in the numerators of the partial 
fraction decomposition. Therefore, the quasipolynomial given by the first 
algorithm has the same value as the quasipolynomial given by the sec- 
ond algorithm over ((^^(^ + C*)) ^ Cij for each j. By Lemma [3. II the two 
quasipolynomials must be equal. □ 



^Our algorithm is of complexity 0{ J ) ^ single chamber where N is the number of 

faces and n is the dimension of the vector space. This is practically applicable (and realized in 
the "vector partition function" program) on a 32-bit machine (low-class modern PC) for more 
than a thousand chambers only if each has up to about 30 faces in S*** dimension, provided 
we run the "vector partition function" program for hours. 
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5 On the periods of the Kostant parti- 
tion function of Eq, Ej, E^, F4 and G2 

The periods for the Kostant partition function for the classical types root 
systems have been computed in [BBC V06] . For type A the period is one 
(i.e. "no period"), for types B, C and D - 2. 

Proposition 5.1 

(a) We are using Definition The period of the Kostant partition 
function of a simple root system A is a divisor of the least com- 
mon multiple of all possible orders of the elements of the group 
Az(A)/Az(Ai), where Ai runs over all root subsystems of A that 
have the same rank as A. 

(b) The possible orders and the upper bounds for the Kostant partition 
function periods are given in the below table. 

Root system 

An Bn Cn Dn E(i Eg F4 G2 

Possible orders of elements in Az(A)/Az(Ai) 

1 1,2 1,2 1,2 1,2,3 1,2,3,4 1,2,3,4,5,6 1,2,3,4 1,^,3 
Kostant function period is a divisor of 

1 2 2 2 6 12 60 12 6 

Proof. 

(a) A priori, the algorithm in Section 12.31 gives us a partial fraction 
decomposition for which the fractions in ^ have fraction support (see 
Section [23] for definition effraction support) consisting not necessarily of 
roots, but of root multiples instead. Ifowever for a root system, accord- 
ing to Lemma [5.21 below, in Step 1, Section [2.31 of our algorithm, we can 
always choose our "preferred linear combination" such that, after even- 
tual re-indexing, ao — 1. Therefore we can always get a partial fraction 
decomposition with denominators corresponding to roots. 

Let ai , . . . a„ be n linearly independent roots, where n is the rank of 
the root system, and let Ai be the root system generated by the Oi's. 

Let n, {t^) = E^e A,(A) ''"ds- 1,0,^(7)2:^) where Az(A) is the lattice 
generated by our root system, and where in the expression Tjs-i_o,d we 
have that B is the matrix (a^), stands for the zero vector and d is 
the corresponding period. By the remarks after ([6]), d is the least com- 
mon multiple of the orders of the elements of Az(A)/Az{Ai}. By ([5| we 
see that the period corresponding to the partial fraction decomposition 
summands involving the roots ai is d. 

(b) Let Ai be a root subsystem of A of the same rank as A. Let 
a £ Az(A)\Az(Ai) be a representative of an element of Az(A)/Az(Ai). 
Modifying a by elements of Ai does not change the represented element, 
therefore we can assume that a is Ai-lowest. All possibilities for Ai- 
lowest vectors for the exceptional Lie algebras are given in the tables of 
Lemma 15.21 below, and by the linear relations there, we see that their 
orders are as in the given table. If A is of type A then Ai = A by 
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the rank Ai = rank A requirement. We leave cases B, C and D to the 
reader with the note that by the rank requirement and [Dyn72[ Table 9] 
Ai cannot have components of type A. □ 

The following lemma is key to Proposition 15.11 It can be derived as 
a consequence of |FT07I Lemma 4]. The present proof does not rely on 
[FT07 and was discovered independently by the author. 

Lemma 5.2 Let A"'" be the positive root system of a simple root system 
of arbitrary type (A,B,C,D,E,F or G) of rank n. Let ai, . . . ,ak € A^ be 
k linearly dependent elements o/ A"*". Then there exists a linear relation 
between the Oi for which one of the roots is with coefficient +1 or -1. 

Proof. Throughout this proof, whenever we are considering exceptional 
Lie algebras (of types Eg, E7, Es, F4 and G2), the coordinates are picked 
in a simple root basis so that the corresponding Cartan matrices are as in 
[Hum72l Table l,page 59], see 

|http: //www. liegroups . org/dissemination/spherical/explorer/rootS ystem. cgi| 

for actual printouts of the root systems. 

Take s linearly dependent vectors, say ai, . . . ,as such that any s — 1 
of them are linearly independent. Take Oi £ Z, such that 

aiQi H + OsOs = (18) 

and so that Os is of minimal possible absolute value (i.e. choose a minimal 
linear relation). Let Ai be the root system generated by ai, . . . , Os-i, and 
let Pi, . . . , /3a-i be a simple basis of Ai with respect to a chosen regular 
height function h, which we will specify only for type A and the exceptional 
Lie algebras. Then the roots Pi are integral linear combinations of the 
roots Oi, . . . , Qs. (|18p becomes 

H \-bs^iPs^i+a,as ^0 (19) 

for some integers bi. Since asQs — as{as — Pi) + asA we see that the 
possible values of a^ do not change if we change ois by the Pi, so we can 
assume that Qs is the lowest weight with respect to Ai (i.e. that Os — Pi is 
not a root Vi). We note that the collection of roots Ai U{as} cannot split 
since otherwise (|18p would not be minimal and let us call A the simple 
root system generated by Ai U {cts}- |Dyn72| Table 9] gives all possible 
maximal (but not necessarily of max;imal rank) root subsystems up to 
conjugation for Ai when A is of type A, B, C and D (note that there are 
two cases for type D) and [Dyn72| Table 10]- when A is of exceptional 
type. We will reproduce the tables within our case by case study. We will 
treat the the case for A of type A separately from types B, C and D, and 
we will treat each exceptional case in its own right. 

• A is of type As-i- Ai is then of type A^-^ + ■ • ■ + Aki with '^^{kr + 
1) = s, and the requirement that Ai be of rank s — 1 forces Ai to 
be of type As-i. Choose a regular height function h so that Qs is 
the /i-lowest weight. Then Qs + /3i + ■ ■ ■ + Ps-i = and we have the 
desired linear combination with — ±1. 

• A is of type Cs-i or Ds-i. Then Ai is of type +• ■ • + 
algebras of types B, C and D with "^^{kr + 1)+ the rest of the ranks 
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— s ~ 1. The requirement that A i be of rank s — 1 forces all ki = 0. 
In order to avoid the tedious exploration of all possible cases (and 
there are quite a few), we will make a Cartan matrix determinant 
trick. Let (•, •)k denote the scalar product on the roots induced by 
the Killing form and normalized so that the long roots have length 
V2 for types B and D and length 2 for type C (this way the most 
predominant root length is -\/2)- For the time being, assume that 
H19p is a linear combination over Q (instead of over Z). To avoid 
confusion rename temporarily the coefficients in (I19p to 2s, bi. Take 
the scalar products {l3i,»)K to obtain 



Now M is a matrix corresponding to the simple basis . . . , /3s-i, 
which is block-diagonal (up to permuting the simple basis vectors 
Pi) with blocks of types B, C and/or D. The determinants of those 
matrices are proportional to the determinants of the Cartan matrices 
depending on the choice (normalization) of {•, For our partic- 
ular choice of normalization, the determinants in question are 1 for 
type B and 4 for types C and D. Therefore, the denominators in the 
rational numbers bi are 1, 2 or 4. Thus we can rescale (|20p to an 
integral linear combination for which as is a divisor of 4. So far we 
have used nothing specific about our choice of the index s; therefore 
we can apply the same reasoning to obtain that ai is a divisor of 4 
for all i. Now at least one ai must be ±1 because otherwise all co- 
efficients in our linear relation H18|) would be even, in contradiction 
with its minimality. 

• A is of type E^. All possibilities for Ai 7^ A come in two conjugacy 
classes ( |Dyn72| Table 10]). Since all roots of Ea are of equal length, 
we can choose the regular height function h on the roots so that ar 
is a Ai-lowest weight (i.e. 0:7 — /? is not a root for any /i-positive root 
P). We present the Dynkin diagrams of the two conjugacy classes 
for Ai 7^ A, an explicit simple coordinate realization of both, a list 
of all possibilities for ar, and an explicit linear relation giving us the 
possible values for ay and the bi. 




Now set as = 1 and solve for the rest of the bi's. Then 
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Case 3^42 






) 


: 


/3i = (0, 0, 0, 0, 0, 


1) 








/32 = (1, 1, 1, 1, 1, 


0) 




) 


: 


/33 = (1, 0, 1, 1, 0, 


0) 








/34 = (0, 1, 1, 2, 2, 


1) 






: 


/35 = (0, 0, 1, 0, 0, 


0) 




) 




/36 = (0, 1, 0, 1, 0, 


0) 




Ai-lowest a^ 


Linear combination 






(-1,-1,-2, 


-2,-1,-1) 


3a7 + 2/3i + /?2 -f 2/?3 + /?4 -f 2/?5 


-f /36 


= 


(-1,-2,-2, 


-3,-2,-1) 


3a7 /?i + 2/?2 + /?3 4- 2/?4 + /?5 + 2/36 


= 






Case ^5 + Ax 








1 


fix= (1, 0, 0, 0, 0, 


0) 






2 


/32= (0, 0, 1, 0, 0, 


0) 






3 


/33= (0, 0, 0, 1, 0, 


0) 






4 


/34= (0, 0, 0, 0, 1, 


0) 






5 


ft = (0, 0, 0, 0, 0, 


1) 




c 


6 


/36= (1, 2, 2, 3, 2, 


1) 




Ai-lowest OL-j 


Linear combination 






(-1,-1,-2,- 


3,-2, -1) 


2ar + /3i + 2/32 + 3/33 2^4 + /35 + /36 


= 



Suppose Ai is of type 3^2- Then qi, . . . ,a6 must each belong to 
exactly one copy of Ai. By the following mini- lemma we get the 
desired linear combination with a coefficient ±1 for some a;. 
Mini-Lemma. Let /3i, /32 be a simple basis of A2. Suppose 01,02 
are roots of A2 such that /3i and /32 are linear combinations of ai and 
ai. (a) Let /Ji -f 2/32 = aioi +a2Q2. (b) Let ±/3i ±/32 = aiQi -1-0202- 
Then either ai or a2 is ±1. 

Proof, (a) Depending on the angle between a\ and Q2 we get 
either /3i = oi or 2 and /32 = 02 resp. 1 or /Ji or 2 = ±(ai — 02) and 
/32 reap. 1 = 02 or 1- lu all possible cases, we get either a\ or ai to be 
±1. We leave (b) to the reader. □ 

If for some ai the root system As-i generated by aj,j 7^ z, is of 
type ?>Ai we are done by repeating the preceding considerations for 
the index i. Similarly we are done if Ag_i is of type E(, - then ai 
is an integral linear combination of the rest of the roots. Suppose 
finally As-i are all of type A5 + Ai. Then choosing difi'erent height 
functions we see from the second table that all ai must be equal to 
2, in contradiction with the minimality of (|18|l . 

• A is of type E7. All possibilities for Ai 7^ A come in six conjugacy 
classes ( |Dyn72| Table 10]). Since all roots of E7 are of equal length, 
we can choose the regular height function h on the roots so that ag 
is a Ai-lowest weight (i.e. ag — /3 is not a root for any /i-positive 
root /3). 
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Case A2 + Ar, 














i: 






/Ji = (l, 0, 


0, 


0, 


0, 


0, 


0) 










02 = {0, 0, 


1, 


0, 


0, 


0, 


0) 




? ^ 

1 4 






/53 = (0, 0, 


0, 


0, 


1, 


0, 


0) 










/34 = (0, 0, 


0, 


0, 


0, 


1, 


0) 




? ^ 






/35 = (0, 0, 


0, 


0, 


0, 


0, 


1) 




? ^ 






^6 = (1, 1, 


2, 


3, 


2, 


1, 


0) 




i 7 






07 = (0, 1, 


0, 


0, 


0, 


0, 


0) 


Ai- 


lowest as 






Linear combination 






(-1,-1, 


-2,-2,-2,-2, 


1) 


Sqs - 


h /3i + 2/32 2/3: 


-h4/: 


4 -f 3As 


f 2/3 


6+07=0 


(-2,-2, 


-3,-4,-3,-2, 


1) 


Sas - 


h 2/3i + /32 /33 


+ 2/34 


+ 


3/?5 + 


4:136 


+ 2/37 = 



By direct observation our mini-lemma from the Ee case applies to 
the vectors 0i and 02 and since two of the vectors a, must be in the 
A2 component we are done. 







Case De + Ai 










/3i = (l, 0, 0, 0, 0, 


0, 


0) 






/32 = (0, 1, 0, 1, 1, 


1, 


1) 




X 


03 = {0, 0, 1, 0, 0, 


0, 


0) 






04 = {0, 0, 0, 1, 0, 


0, 


0) 






/3b = (0, 0, 0, 0, 1, 


0, 


0) 






06 = {0, 0, 0, 0, 0, 


1, 


0) 




7 


07 = {1, 2, 2, 3, 2, 


1, 


0) 


Ai- 


lowest as 


Linear combination 






(-2,-2, 


■3,-4,-3,-2,-1) 


2as + 3/3i -1- 2/32 + 4/33 + 304 + 2^b 


+ 


6+07 = 



Case ^3-1-^1-1- A3 



Ai-lowest as 
(-1,-1,-2,-3,-2,-2,-1) 
(-2,-2,-3,-4,-3,-2,-1) 
(-1,-1,-2,-2,-1,-1,-1) 



01 


= (1, 


0, 


0, 


0, 


0, 


0, 


0) 


02 


= (0, 


0, 


1, 


0, 


0, 


0, 


0) 


/33 


= (0, 


0, 


0, 


1, 


0, 


0, 


0) 


/34 


= (1, 


1, 


2, 


3, 


3, 


2, 


1) 


/35 


= (0, 


0, 


0, 


0, 


0, 


1, 


0) 


06 


= (0, 


0, 


0, 


0, 


0, 


0, 


1) 


07 


= (1, 


2, 


2, 


3, 


2, 


1, 


0) 



Linear combination 
4a8 +I3i + 2/32 + 3/33 + 2/34 + Ws, 2/36 + /^t = 
4as + m + 202 + P?, + 2/34 + 05 + 206 + 3/37 = 
2as + 01+ 2/32 4- /33 -f /35 -f 2/36 + /37 = 
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Case A7 



Pi 


= (1, 


0, 


0, 


0, 


0, 


0, 


0) 


P2 


= (0, 


0, 


1, 


0, 


0, 


0, 


0) 


/33 


= (0, 


0, 


0, 


0, 


1, 


0, 


0) 


Pa 


= (0, 


0, 


0, 


0, 


0, 


1, 


0) 


P5 


= (0, 


0, 


0, 


0, 


0, 


0, 


1) 


Pe 


= (1, 


1, 


2, 


3, 


2, 


1, 


0) 


Pr 


= (0, 


1, 


0, 


0, 


0, 


0, 


0) 



Ai-lowest as 
(-1,-1,-2,-3,-3,-2,-1) 



2as 



Linear combination 
- Pi + 2/32 + 3^3 + 4^4 + SPs + 2P6 + P7 = 



Case D4, + 3Ai 




Pi 


= (1, 


0, 


0, 


0, 


0, 


0, 


0) 


P2 


= (0, 


0, 


0, 


1, 


0, 


0, 


0) 


P3 


= (0, 


0, 


1, 


0, 


0, 


0, 


0) 


Pi 


= (0, 


1, 


0, 


1, 


1, 


1, 


1) 


Ps 


= (1, 


1, 


2, 


3, 


3, 


2, 


1) 


Pe 


= (0, 


0, 


0, 


0, 


0, 


1, 


0) 


P7 


= (1, 


2, 


2, 


3, 


2, 


1, 


0) 



Ai-lowest Qg 
(-2,-2,-3,-4,-3,-2,-1) 
(-1,-1,-2,-3,-2,-2,-1) 
(-1,-2,-2,-3,-2,-2,-1) 



Linear combination 

2Qg + 2/3i +132+ 2/33 +Pi+P5+P7 

2q8 + Pi + 2/32 + 2/33 + P4 + P5+P6 
2as + P1 + P2 + 2/33 + 2/34 + Pe + P7 





Case 7Ai 












1 

2 


Pi = {l, 


0, 


0, 0, 


0, 




0) 




P2 = {0, 


0, 


1, 0, 


0, 




0) 


o3 


P3 = (0, 1 


0, 


1, 1, 


1, 




1) 


4 


Pi = (1, 1 


2, 


3, 3, 


2, 




1) 


5 


/35 = (0, 


0, 


0, 0, 


1, 




0) 


6 


Pe = {l, 2 


2, 


3, 2, 


1, 




0) 


7 


P7 = (1, 1 


2, 


2, 1, 


1, 




1) 


Ai-lowest Qg 


Linear combination 








(-1,-1,-1,-2,-1,-1,0) 


2Qg + Pi 


+ P2 


+ Pr. + Pa 









(-1,-1,-1,-2,-1,-1,-1) 


2Qg + Pi 


+ P2 


+ P3 + P7 









(-1,-1,-1,-2,-2,-2,-1) 


2Qg + Pi 


+ P3 


+ Pi + P5 









(-1,-1,-2,-3,-2,-2,-1) 


2as + P2 


+ Pi 


+ P5 + P7 









(-1,-2,-2,-3,-2,-2,-1) 


2Qg + Pa 


+ P5 


+ P6 + P7 









(-1,-2,-2,-4,-3,-2,-1) 


2a8 + P2 


+ P3 


+ Pi + Pa 









(-2,-2,-3,-4,-3,-2,-1) 


2Qg + Pi 


+ Pi 


+ P6 + P7 










If for some q; the root system Ag-i generated by aj,j 7^ i, is of type 
A2 + Ac, we are done by repeating the considerations from the very 
first case. Similarly we are done if Ag-i is of type E7 - then Qi is an 
integral linear combination of the rest of the roots. Suppose finally 
Ag-i are all not of type A5 + A2 or £7 and assume Oi is Ag^^-lowest. 
Then (possibly choosing different height functions) we see from the 
remaining 5 tables that all at must be even, in contradiction with 
the minimality of (|18p . 
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• A is of type -Eg- All possibilities for Ai ^ A come in fourteen 
conjugacy classes ( [Dyn72[ Table 10]). Since all roots of -Eg are of 
equal length, we can choose the regular height function h on the 
roots so that og is a Ai-lowest weight (i.e. og — /3 is not a root for 
any /i-positive root /3). 



Case As 



/3i 


= (1, 


3, 


3, 


5, 


4, 


3, 


2, 


1) 


/32 


= (1, 


0, 


0, 


0, 


0, 


0, 


0, 


0) 


/33 


= (0, 


0, 


1, 


0, 


0, 


0, 


0, 


0) 


/34 


= (0, 


0, 


0, 


1, 


0, 


0, 


0, 


0) 


/35 


= (0, 


0, 


0, 


0, 


1, 


0, 


0, 


0) 


/36 


= (0, 


0, 


0, 


0, 


0, 


1, 


0, 


0) 


/37 


= (0, 


0, 


0, 


0, 


0, 


0, 


1, 


0) 


/38 


= (0, 


0, 


0, 


0, 


0, 


0, 


0, 


1) 



Ai-lowest Qg Linear combination 

(-1,-1,-2,-3,-3,-3,-2,-1) 3a9 4- /3i -h 2/32 + + 4/^4 + 5/35 + 6/36 4/37 + 2/3g = 
(-2,-2,-4,-5,-4,-3,-2,-1) 3a9 2/3i + 4/32 + 6/33 5/34 -t- 4/35 + 3/36 2/37 -f /3g = 



Ag has 36 positive roots and there are only 




= 30260340 



possible choices of roots ai, . . . Qg. Direct computation of all possi- 
bilities shows that for all such choices of maximal rank, there is at 
least one on with coefficient ±]|f|. 





Case 2A4 














9 ^ 
1 2 


/3i = (0, 0, 


0, 


0, 


1, 0, 


0, 


0) 






/32 = (0, 0, 


0, 


1, 


0, 0, 


0, 


0) 




? ^ 


/33 = (0, 0, 


1, 


0, 


0, 0, 


0, 


0) 






/34 = (1, 0, 


0, 


0, 


0, 0, 


0, 


0) 




1 ^ 


/35 = (1, 1, 


2, 


3, 


3, 3, 


2, 


1) 




if 6 


/36 = (1, 2, 


2, 


3, 


2, 1, 


0, 


0) 




1 r 


/37 = (0, 0, 


0, 


0, 


0, 0, 


1, 


0) 




I 8 


/38 = (0, 0, 


0, 


0, 


0, 0, 


0, 


1) 




Ai-lowest Qg 


Linear combination 








(-1,-1,-2,-3,-2,-1,0,0) 


5q9 + 3/3i 6/32 + 4/3: 


+ 2/34 /35 2^6 


- 2/37 - Ps 


= 


(-1,-1,-2,-3,-2,-1,-1,-1) 


5Qg 3/3i + 6/32 4/33 


+ 213 


i + l35 + 2/36 


+ 3/37 + 4/3g 


= 


(-1,-1,-2,-3,-3,-2,-1,0) 


5a9 + 4^1 + 3/32 4- 2^ 


3+/3 


4 4- 3/35 + /36 


-/37 


-3/3g 


= 


(-1,-1,-2,-3,-3,-2,-1,-1) 


5ag -f 4^1 -1- 3/32 + 2^ 


3+/3 


4 3/35 + /36 




+ 2/3g 


= 


(-2,-2,-3,-4,-3,-2,-1,-1) 


5q9 -f /3i -f 2/32 + 3/33 


+ 4/34 


-h 2/35 + 4/36 


+ 13 


7+3/38 


= 


(-2,-2,-3,-4,-3,-2,-2,-1) 


5q9 /3i 2/32 + 3/33 + 4/34 


+ 2/35 + 4/36 


+ 6/37 + 3/3g 


= 


(-2,-2,-4,-5,-4,-3,-2,-1) 


5q9 + 2Pi + 4/32 + 6/33 


+ 3/3 


4 4/35 3/3 


6 2/37 /3g 


= 


(-2,-3,-4,-6,-5,-3,-2,-1) 


5a9 4^91 + 3/32 2/33 


+ /34 


3/35 + 6/36 


+ 4/; 


(7 + 2/3g 


= 



^ All easy considerations, attempted by tlie author and similar to the ones given for the 
i?7 and Es cases, failed in the current case. We therefore enumerated all possible 30260340 
combinations by brute computer force and verified the claim (the computation itself took less 
than 45 minutes). 
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2A4 has 20 positive roots and there are only 



10 

4 



: 44100 possi- 



ble choices of roots qi, . . . ag. Direct computation of all possibilities 
shows that for all such choices of maximal rank, there is at least one 
Qi with coefBcient ifl 



Case A5 + Ai+ A2 



Pi 


= (0, 


0, 


0, 


0, 


0, 


1, 


0, 


0) 




= (0, 


0, 


0, 


0, 


1, 


0, 


0, 


0) 


/33 


= (0, 


0, 


0, 


1, 


0, 


0, 


0, 


0) 


/34 


= (0, 


0, 


1, 


0, 


0, 


0, 


0, 


0) 


/35 


= (1, 


0, 


0, 


0, 


0, 


0, 


0, 


0) 




= (1, 


2, 


2, 


3, 


2, 


1, 


0, 


0) 




= (2, 


3, 


4, 


6, 


5, 


4, 


3, 


1) 


Ps 


= (0, 


0, 


0, 


0, 


0, 


0, 


0, 


1) 



Ai-lowest Qg 
(-1,-1,-2,-3,-2,-1,0,0) 
(-1,-1,-2,-3,-3,-2,-1,0) 
(-1,-1,-2,-3,-3,-2,-1,-1) 
(-2,-2,-3,-4,-3,-2,-1,-1) 
(-2,-2,-4,-5,-4,-3,-2,-1) 
(-2,-3,-4,-6,-5,-4,-2,-1) 



Linear combination 
209 + Pi + 2/?2 + 3/33 + 2Pi + Ps + Pa = 
3a9 + 2/?l + 4/?2 + 3/?3 + 2/?4 + /?5 + /?7 - /?8 = 
Sag + 2/3i + 4/32 + 3/33 + 2/34 + P5 + Pr + 2Ps ^ 
6ag + 2/32 + 3/33 + 4/34 + 5/35 + 3/36 + 2/37 + 4/38 = 

3ag + ^1 + 2/32 + 3/33 + 4/34 + 2/35 + 2/37 4- /38 = 
Gag + 5/3i 4/32 + 3/33 2/34 + /35 + 3/3^ + 4/37 + 2/38 = 



For the cases for which the coefficient in front of ag is 3, the mini- 
lemma from the Ea case applies to the roots ^7 and Ps to produce 
a relation in the ai's with coefficient ±1. The case when the coeffi- 
cient in front of ag is even is treated together with the most of the 
remaining cases around at the end of the Eg proof. 



^We have, for a second time, summoned our computer computational power. 
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Case A2 + Ed 



I3i 


= (2, 


3, 


4, 


6, 


5, 


4, 


3, 


1) 


P2 


= (0, 


0, 


0, 


0, 


0, 


0, 


0, 


1) 




= (1, 


0, 


0, 


0, 


0, 


0, 


0, 


0) 




= (0, 


1, 


0, 


0, 


0, 


0, 


0, 


0) 


/35 


= (0, 


0, 


1, 


0, 


0, 


0, 


0, 


0) 


/36 


= (0, 


0, 


0, 


1, 


0, 


0, 


0, 


0) 




= (0, 


0, 


0, 


0, 


1, 


0, 


0, 


0) 


/Js 


= (0, 


0, 


0, 


0, 


0, 


1, 


0, 


0) 



Ai-lowest ag Linear combination 

(-1,-1,-2,-3,-2,-1,-1,-1) 3a9 + + 2/32 + Pi + 2/35 3/36 -I- /37 - /38 = 

(-1,-1,-2,-2,-2,-2,-1,-1) Sag + Pi + 2/32 + + 2/35 + P7 + 2Ps=0 

(-1,-1,-2,-3,-3,-2,-1,0) Sag + Pi - P2 + + 2p5 + SPe + 4/37 + 2Ps ^ 

(-1,-1,-2,-3,-3,-2,-1,-1) 3^9 + Pi+ 2/32 + Ps + 2/35 + SPe + 4/37 + 2/3^ = 

(-1,-2,-2,-3,-3,-2,-1,-1) 3a9 + Pi + 2/32 + Pa + 3/34 + 2/35 + 3/36 + 4/37 + 2/3s = 

(-2,-2,-3,-4,-3,-2,-1,-1) 3a9 + Pi + 2P2 + 4/33 + ?>Pa + 5/35 + 6/36 -f 4/37 4- 2/38 = 

(-2,-2,-3,-4,-3,-2,-2,-1) Sag + 2/3i + P2 + 2^3 + Pr> - Pi - 2Ps = Q 

(-2,-2,-3,-4,-4,-3,-2,-1) Sag + 2Pi + P2 + 2Pz + Pz + 2Py + Ps = Q 

(-2,-2,-4,-5,-4,-3,-2,-1) 3^9 -I- 2/3i + P2 + 2p3 + 4/35 + Sp6 + 2p7 + ps = 

(-2,-3,-4,-5,-4,-3,-2,-1) 3q9 + 2/3i + P2 + 2/33 + + 4^5 + 3/36 + 2/37 + /3s = 

(-2,-3,-4,-6,-5,-3,-2,-1) 309 + 2/3i + P2 + 2/33 -H 3/34 + 4/35 -H 6/36 -h 5/37 -h /3s = 

(-2,-3,-4,-6,-5,-4,-2,-1) 3a9 + 2,/3i + & + 2/33 + 3/34 + 4/35 + 6/36 H- 5/37 4/38 = 



The mini-lemma from the Eq case applies to the roots Pi and /32 to 
produce a relation that contains a coefficient ±1. 
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Case 4^2 














/9i = (0, 0, 


0, 


0, 


0, 


1, 


0, 


0) 


02 = {0, 0, 


0, 


0, 


1, 


0, 


0, 


0) 


03 = (1, 1, 


2, 


3, 


2, 


1, 


0, 


0) 


04 = (0, 1, 


0, 


0, 


0, 


0, 


0, 


0) 


05 = {0, 0, 


1, 


0, 


0, 


0, 


0, 


0) 


06 = {1, 0, 


0, 


0, 


0, 


0, 


0, 


0) 


^7 = (2, 3, 


4, 


6, 


5, 


4, 


3, 


1) 


^8 = (0, 0, 


0, 


0, 


0, 


0, 


0, 


1) 



Ai-lowest ag Linear combination 

(-1,-1,-1,-1,-1,-1,0,0) Sag -I- 2/3i +^2+03 + 2/34. +05 + 2/36 = 

(-1,-1,-2,-2,-2,-1,0,0) 3a9 + 0i + 202 + 2/33 + 04 + 2/3^+06=0 

(-1,-1,-2,-3,-2,-1,-1,-1) 3q9 - 2l3i -132+ 0a - /34 4- /37 -f 2/38 = 

(-1,-1,-2,-2,-2,-2,-1,-1) 3q9 + 2/3i +132 + 2/3,, + /36 + 07 + 2/3^ = 

(-1,-1,-2,-3,-3,-2,-1,0) 3a9 -f /3i 2/32 -f /^s - /34 /37 - /^s = 

(-1,-1,-2,-3,-3,-2,-1,-1) 3a9 + 0i + 202 -h /Ss - ^4 /37 2/38 = 

(-1,-2,-2,-3,-3,-2,-1,-1) 3a9 + 0i + 202 + 03 + 204. +07 + 2/38 = 

(-2,-2,-3,-4,-3,-2,-1,-1) 3q9 + 2/3., +04 + 05 + 206 +07 + 2/3s = O 

(-2,-2,-3,-4,-3,-2,-2,-1) 3q9 - 2/3i - /32 + 05 + 2/36 + 2P7+/3s = 

(-2,-2,-3,-4,-4,-3,-2,-1) 3a9 + 0i + 202 + 05 + 2/36 +207+0^ = 

(-2,-2,-4,-5,-4,-3,-2,-1) Sag +03 - 04 + 2/35 + 06 + 2/37 + I3s = 

(-2,-3,-4,-5,-4,-3,-2,-1) Sag +03 + 204. + 205 -h /Js -h 2^7 -h /Ss = 

(-2,-3,-4,-6,-5,-3,-2,-1) 3q9 - 0i + 02 + 2/33 /34 2/37 /38 = 

(-2,-3,-4,-6,-5,-4,-2,-1) 3a9 + 2/3i + /32 + 2/33 + 04 + 207 + 0ii = 



The mini-lemma from the Ee case apphes to the roots /3i and 02 to 
produce a relation that contains a coefficient ±1. 







Case Ds 














1 


o 2 


01 = (1, 


2, 


2, 


3, 


2, 


1, 


0, 


0) 




/3 


02 = (0, 


0, 


0, 


0, 


0, 


0, 


0, 


1) 




4 


03 = (0, 


0, 


0, 


0, 


0, 


0, 


1, 


0) 




5 


04 = (0, 


0, 


0, 


0, 


0, 


1, 


0, 


0) 




6 


05 = (0, 


0, 


0, 


0, 


1, 


0, 


0, 


0) 




7 


06 = (0, 


0, 


0, 


1, 


0, 


0, 


0, 


0) 




8 


07 = (0, 


0, 


1, 


0, 


0, 


0, 


0, 


0) 






0i = (1, 


0, 


0, 


0, 


0, 


0, 


0, 


0) 


Ai-lowest ag 




Linear combination 






(-2,-3,-4,-6, 


-5,-4,-3,-2) 


2ag + 301 + 402 + 603 


+ 50 


4 + 4:05 


+ 30 


3 -f 2/37 -h /38 = 
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Case Ai + Ai 



o 1 


/3i = (2, 


3, 


4, 


6, 


5, 


4, 


3, 


2) 


9 ^ 


01 = (1, 


2, 


2, 


3, 


2, 


1, 


0, 


0) 


I 3 


h = (0, 


0, 


0, 


0, 


0, 


0, 


1, 


0) 


I 4 


/34 = (0, 


0, 


0, 


0, 


0, 


1, 


0, 


0) 


I 5 


/35 = (0, 


0, 


0, 


0, 


1, 


0, 


0, 


0) 


I 6 


/36 = (0, 


0, 


0, 


1, 


0, 


0, 


0, 


0) 


I 7 


^7 = (0, 


0, 


1, 


0, 


0, 


0, 


0, 


0) 


I 8 


08 = (1, 


0, 


0, 


0, 


0, 


0, 


0, 


0) 



Ai-lowest at) Linear combination 

(-1,-1,-2,-3,-3,-2,-1,0) 2a9 + 02 + 2(33 + i0A + 4:^5 + S06 + 2^7 + 0s = 

(-2,-2,-4,-5,-4,-3,-2,-1) 4q9 + 2/3i + /32 + 2^3 + 304 + 4^b + 506 + 607 + 30a = 

(-2,-3,-4,-6,-5,-4,-3,-1) 4q9 + 2/3i + 302 + 603 + 504 + 4/35 + 306 + 2/37 + 08 =0 







Case Ai + E7 














c 


1 


01 = {2, 3, 


4, 


6, 


5, 


4, 


3, 


1) 






2 


02 = {1, 0, 


0, 


0, 


0, 


0, 


0, 


0) 






4 


/?3 = (0, 1, 


0, 


0, 


0, 


0, 


0, 


0) 






-S_o 3 


/34 = (0, 0, 


1, 


0, 


0, 


0, 


0, 


0) 






6 


05 = {0, 0, 


0, 


1, 


0, 


0, 


0, 


0) 






7 


06 = {0, 0, 


0, 


0, 


1, 


0, 


0, 


0) 






8 


07 = {0, 0, 


0, 


0, 


0, 


1, 


0, 


0) 








^8 = (0, 0, 


0, 


0, 


0, 


0, 


1, 


1) 




Ai-lowcst ag 


Linear combination 








(-1,-1,-2,-3 


-3,-2,-1,0) 


2ag+l3 


1-03 


+ / 


36- 


'08 = 


-- 






(-2,-2,-3,-4, 


-3,-2,-2,-1) 


2a9 +I3i + 2/?2 + 03 + 2/34 


+ 2A5 + 0<i 


+ 08=0 




(-2,-2,-3,-4, 


-4,-3,-2,-1) 


2a» + I3i + 2/32 + 03 


+ 2/34 


+ 2/35 + 3/36 


+ 2/37 + 08 


= 


(-2,-2,-4,-5, 


-4,-3,-2,-1) 


2aa + Bi + 2/32 + 03 + 4/34 


+ 4/3,5 + 3/36 


+ 2/37 + 08 


= 


(-2,-3,-4,-5, 


-4,-3,-2,-1) 


2q9 + 01+ 2/32 + 3/33 


+ 4/34 


-H4/35 


+ 3l3e 


+ 2/37 + 08 


= 


(-2,-3,-4,-6,-5,-3,-2,-1) 


209 +01 + 202 + 303 


+ 4^4 


-H6/35 


+ 50e 


+ 207 + 08 


= 


(-2,-3,-4,-6,-5,-4,-2,-1) 


2ag + (3i+ 2/32 + 303 


+ 4134 


+ 6/35 -f 5/3c 


+ 4/37 + 08 


= 


(-2,-3,-4,-6, 


-5,-4,-3,-2) 


2q9 +01 + 2/32 + 303 


+ 4/34 


6/35 5/36 


+ 4/37 + 3/38 


= 
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Case De + 2Ai 










ol o2 


01 = (1, 


0, 0, 


0, 0, 0, 


0, 


0) 






02 = (0, 


1, 0, 


1, 1, 1, 


1, 


0) 






03 = (0, 


0, 1, 


0, 0, 0, 


0, 


0) 






0i = (0, 


0, 0, 


1, 0, 0, 


0, 


0) 




a 


05 = (0, 


0, 0, 


0, 1, 0, 


0, 


0) 




o 7 


06 = (0, 


0, 0, 


0, 0, 1, 


0, 


0) 




o 8 


07 = (2, 


3, 4, 


6, 5, 4, 


3, 


2) 






/38 = (1, 


2, 2, 


3, 2, 1, 


0, 


0) 




^l-loW6St CK9 


Linear combination 








(_1._1,_1._1,_1._1,0,0) 


2a9 + /3i - /34 


+ 06+08- 


- 






f_l_l -2-2,-2,-1,0,0) 


2a9 /3i 2/33 + /34 


+ 205 + 06 


+ 0i 


= 




f-1 -1 -2 -3 -2 -1 0) 


2q9 -h /3i 2/33 3/3 


1+2^5+ 06 


+ 


8=0 




(-1^-1.-2,-3,-2,-1,-1,-1) 


2a9 - 


-02+04 


- 06 + 07 - 


= 






(-1,-1,-2,-2,-2,-2,-1,-1) 


2a9 - 


-02-04 


+ 06 +07 - 


= 






(-1,-1,-2,-3,-3,-2,-1,-1) 


2q9 - /3: 


+ I3i+ 2/35 + /36 + /37 = 







(-1,-2,-2,-3,-3,-2,-1,-1) 


2a9 -01- 


02 - 2/33 


- 2/34 -1- 07 


+0i 


i = 




(-2,-2,-3,-4,-3,-2,-1,0) 


2q9 3/3i + 2/32 


-1-4^3-1- 3^84-1-2/35-1-/36-1-/38 = 


= 


(-2,-2,-3,-4,-3,-2,-1,-1) 


2a9 + /3i - 02 


+ 07 + 08 - 


- 






(-2,-2,-3,-4,-3,-2,-2,-1) 


2af) + 2/3i + 132 + 2/33 + /34 - /36 


+ 0- 


' =0 




(-2,-2,-3,-4,-4,-3,-2,-1) 


2a9 + 2/3i -h /3: 


+ 2/33 + /34 + 2,/35 + 


06 H 


-07 = 





(-2,-2,-4,-5,-4,-3,-2,-1) 


2a9 + 2Pi + 132 


+ 4/33 + 3/34 -1- 2/35 -1- /36 


+ 07 = 





(-2,-3,-4,-5,-4,-3,-2,-1) 


2a9 -1- 01 


-1- ft -1- 2/33 +07 + 08 = 







(-2,-3,-4,-6,-5,-3,-2,-1) 


2ag +01+^2 


+ 2/33 + 


2/34 + 2/3,5 +07 + 08 = 





(-2,-3,-4,-6,-5,-4,-2,-1) 


2Q9 +01+02+ 2P3 + 2/34 


+ 2/35 + 2/3 


i+07+0i 


! =0 


(-2,-3,-4,-6,-5,-4,-3,-1) 


2q9 + 2/3i -h 3/3: 


-1-4/33 -1-3/34 -1- 2/35 -I-/36 


+ 07 -- 


= 
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Case Ds + A3 




Ai-lowest 09 
-1,-1,-1,-1,-1,-1,0,0) 
-1,-1,-2,-2,-2,-1,0,0) 
,-1,-1,-2,-3,-2,-1,0,0) 
-1,-1,-2,-3,-2,-1,-1,-1) 
-1,-1,-2,-2,-2,-2,-1,-1) 
,-1,-1,-2,-3,-3,-2,-1,0) 
-1,-1,-2,-3,-3,-2,-1,-1) 
-1,-2,-2,-3,-3,-2,-1,-1) 
-2,-2,-3,-4,-3,-2,-1,0) 
-2,-2,-3,-4,-3,-2,-1,-1) 
-2,-2,-3,-4,-3,-2,-2,-1) 
-2,-2,-3,-4,-4,-3,-2,-1) 
-2,-2,-4,-5,-4,-3,-2,-1) 
-2,-3,-4,-5,-4,-3,-2,-1) 
-2,-3,-4,-6,-5,-3,-2,-1) 
-2,-3,-4,-6,-5,-4,-3,-1) 



Pi 


= (1, 


0, 


0, 


0, 


0, 


0, 


0, 


0) 




= (0, 


1, 


0, 


1, 


1, 


1, 


1, 


0) 




= (0, 


0, 


1, 


0, 


0, 


0, 


0, 


0) 




= (0, 


0, 


0, 


1, 


0, 


0, 


0, 


0) 


/35 


= (0, 


0, 


0, 


0, 


1, 


0, 


0, 


0) 


/36 


= (1, 


1, 


2, 


3, 


3, 


3, 


2, 


1) 


07 


= (1, 


2, 


2, 


3, 


2, 


1, 


0, 


0) 


08 


= (0, 


0, 


0, 


0, 


0, 


0, 


1, 


1) 



Linear combination 

■ 2/33 - 4/34 - 2/35 + fie 



4a,j + /3i - /32 - 2/33 - 4/34 - 2/35 + fie + 2/37 - 08 = 

4a9 + + 2/33 + 2/35 + 0e + 207 - 0s = Q 

4Qg + /3i - /32 + 2/33 + 4/34 + 2 fir, + /36 + 2/37 - /38 = 
4a9 + 01-02 + 2/33 + 4/34 + 2/35 + 0(. + 2/37 + 3/38 = 
4a9 - /3i - 3/32 - 2/33 - 4/34 - 2/35 + ifie + 2fi7 + fi»^Q 
2ag + fii+fi2 + 2/33 + 2/34 2/35 /36 - /38 = 
4a9 -fii- 3fi2 - 2/33 + 2/35 + 3/36 2/37 /38 = 
2a9 -01-02- 2fiz - 2fii + 0e + 2fiY + 0% = Q 
4q9 -I- 5/3i -h 3/32 + 6/33 4/34 + 2 fir + fie + 2fi-! - fis = Q 

209 -f /3i - /32 -f /36 -f 2/37 + /38 = 
4q9 + 5/3i 3/32 6/33 4/3* 2/35 + 0e + 2 fir + 30s = 

4a9 + 3/3i /32 2/33 + 2/35 + 306 + 2fi7 + fis = 
4a9 + 301 +fi2+ 603 + 4:04. + 2 fir, + 306 + 2fir + fi8=0 
2ag +fii+fi2 + 2/33 -f /36 -f 2/37 -f /38 = 
2q9 +fii+fi2+ 2fi3 + 2/34 + 2 fir, +fie + 2fi7+fi8 = 
4a9 + 3fii + 502 + 603 + 4/34 + 2/35 + 306 + 2fir+fis=0 
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Case 2Z)4 














/3i = (1, 


0, 


0, 0, 0, 0, 





0) 








/32 = (0, 


1, 


0, 1, 1, 1, 


1, 


0) 




Y 




/33 = (0, 


0, 


1, 0, 0, 0, 


0, 


0) 








/34 = (0, 


0, 


0, 1, 0, 0, 


0, 






o 5 o 6 




/35 = (1, 


1, 


2, 3, 3, 2, 


1, 


0) 




Is 




/36 = (0, 


0, 


0, 0, 0, 1, 


0, 


0) 








= (0, 


0, 


0, 0, 0, 0, 


1, 


1) 








/38 = (1, 


2, 


2, 3, 2, 1, 


0, 


0) 




Ai-lowest cxg 






Linear combination 








(_1^_1^_1^_1^_1^_1^0 


^ } 


2a9-l-/3 


1 - /34 /36 /38 = 


= 






(-1,-1,-2,-2,-2,-1,0 


^ } 


2^9 




2 - /34 /35 -h /38 = 


= 






(-1,-1,-2,-3,-2,-1,0 


^ } 


2a9 




2 /34 /35 -h /38 = 


= 






(_l,_l._2,-3,-2,-l,-l 


-1) 


2Q9 - /32 /34 /35 2^7 /38 = 







(-1,-1,-2,-2,-2,-2,-1 


-1) 


2a9 - /32 - 


-/34 


/35 2/36 -H 2^37 


-H/38 


= 




(-1,-1,-2,-3,-3,-2,-1 


-1) 


2a9 - ,5i - 2/?2 - 


- 2/33 - /34 + 2„55 + i3fi 


+ 2,3 


7 + 38 = 





(-1,-2,-2,-3,-3,-2,-1 


-1) 


2a9 - /3i - /?2 - 


2/33 


- 2/34 + /35 + /36 + 2/37 


+ 2/3s = 





(-2,-2,-3,-4,-3,-2,-1 


,0) 


2a9 2/3i 


+ /? 


2 2/33 /34 /35 


+ 


= 




(-2,-2,-3,-4,-3,-2,-1 


-1) 


2q;9 + /3i - 




+ ,/35 /36 2;37 -H 2/38 


= 




(-2,-2,-3,-4,-3,-2,-2 


-1) 


2a9 2/3i -I- /32 + 


2/33 /34 /35 2/37 /38 = 




(-2,-2,-3,-4,-3,-3,-2 


-1) 


2a9 2;9i /32 2/3 


3 /34 ,35 2^6 


+ 2/3 


7-^,38 = 





(-2,-2,-3,-4,-4,-3,-2 


-1) 


2a9 /3i - 


-/94 


2^5 -h /Se -h 2/37 


+ 138 


= 




(-2,-2,-4,-5,-4,-3,-2 


-1) 


2q9 -I- /3i 2/33 + ^4 -1- 2^B -1- ^6 + 2^7 + /38 = 




(-2,-3,-4,-5,-4,-3,-2 


-1) 


2^9 -h /3i -h /3; 


! 2/33 + /35 /36 2/37 2/38 = 




(-2,-3,-4,-6,-4,-3,-2 


-1) 


2q9 /9i /92 2/33 


2/34 ,35 -h /Se -h 2,37 


-h2,38 = 





(-2,-3,-4,-6,-5,-4,-3 


-1) 


2a9 -h /3i -h 2/32 -h 2/3 


3 -h ,04 -h 2^5 + /36 


-1-2/3 


7+138 = 
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Case D4 + 4Ai 








ol o2 


/3i = (0, 1, 0, 


1, 1, 1, 


1, 


0) 




/32 = (1, 1, 2, 


3, 3, 2, 


1, 


0) 




/33 = (0, 0, 0, 


0, 0, 0, 


0, 


1) 


o 5 


/34 = (1, 1, 2, 


2, 1, 1, 


1, 


0) 


o 6 


/35 = (0, 0, 0, 


1, 0, 0, 


0, 


0) 


o 7 


/36 = (1, 0, 0, 


0. 0, 0, 


0, 


0) 


o 8 


/37 = (1, 2, 2, 


3. 2, 1. 


0. 


0) 




P8 = {0, 0, 0, 


0, 0, 1, 


0, 


0) 


Ai-lowest cxg 


Linear combination 






(_1._1,_1._1,_1._1,0,0) 


2a9 - /35 ,36 


+ 07 + 08 -- 


= 




-2-1 -1,0,0) 


2a9 /35 /36 


+ 07+0S-' 


= 




(-1 -1 _2 -2 -2 -1 01 


2a9 -Pi +132 


-05 +07 -- 


= 




(-1,-1,-2,-3,-2,-1,0,0) 


2ag -pi +132 


+ P5+P7-- 


= 




(-1,-1,-2,-3,-2,-1,-1,-1) 


2a9 ,32 2/33 /34 -h /35 - /38 = 





(-1 1 2,-2,-2,-2,-1,-1) 

V -^5 -^5 -^5 ^1 -^5 ^1 -^5 


2^9 + A + 2/33 +P4-P5+Ps = 





(-1 -2 -2 -3 -3 -2 -1 -1) 


2a9 +I3i + P2 + 2/33 - /35 /37 = 





(-2,-2,-3,-4,-3,-2,-1,0) 

V 7 3 7 3 7 3 7 / 


2a9 + P2+/34. 


+ 06+07 = 


= 




(-2,-2,-3,-4,-3,-2,-1,-1) 


2an +/32 + 2/33 +P4+P6 + P7 = 





(-1,-2,-2,-4,-3,-3,-2,-1) 


2q9 + 2/?i + /32 + 2/3 


3+04+ 05 


+ 0S 


= 


(-2,-2,-3,-4,-3,-2,-2,-1) 


2a9 +P1+P2+ 2/33 


+ 204 + 06 


-0s 


= 


(-2,-2,-3,-4,-3,-3,-2,-1) 


2a9 ,3i ,32 2/33 


+ 2/34 + 06 


+ 0S 


= 


(-2,-2,-3,-4,-4,-3,-2,-1) 


2a9 -h ,3i -1- 2/32 2/3 


3+0i-05 


+ 06 


= 


(-2,-2,-3,-5,-4,-3,-2,-1) 


2q9 + /3i + 2/32 + 2/3 


3 +04 + 05 + 06 


= 


(-2,-3,-3,-5,-4,-3,-2,-1) 


2a9 + 2,/3i + 02 + 2/33 + /34 + /36 


+ 07 


= 


(-2,-3,-4,-5,-4,-3,-2,-1) 


2a9 + 01 + 02 + 2„33 + 2„34 - 0s 


+ 07 


= 


(-2,-3,-4,-6,-4,-3,-2,-1) 


2a9 +P1+P2+ 2/33 


+ 204 + 05 


+ 07 


= 


(-2,-3,-4,-6,-5,-3,-2,-1) 


2q9 +I3i+ 2/32 + 2/3 


3+04+07 


-08 


= 


(-2,-3,-4,-6,-5,-4,-2,-1) 


2q9 -h ,81 -I- 2/32 2,3 


3+04+07 


+ 08 


= 



32 







Case 2A:i + 2Ax 










Pl — (1, 2, 2, 3, 2, 2, 2, 






Is 




/3 /l 1 00001 

Pi = (1, 1, 2, 3, 3, 2, 1, 


0) 








/53 = (0, 0, 0, 0, 0, 0, 0, 


1) 




o 4 

I 




/34 = (0, 0, 0, 1, 0, 0, 0, 


0) 




6 o 




^5 = (0, 0, 1, 0, 0, 0, 0, 
^6 = (1, 0, 0, 0, 0, 0, 0, 


0) 
0) 




o 7 
o 8 




^7 = (1, 2, 2, 3, 2, 1, 0, 
ft = (0, 0, 0, 0, 0, 1, 0, 


0) 
0) 




Al-lowCSt Q9 




Linear combination 






(-1 _1 _1 -1 -1 _1 




2a9 - /34 ft /St + ft = 






f-1 -1 -1 -2 -1 -1 


0) 


2a9 -1- ft -1- ft -1- ft -h ft = 






(_l^_l^_2,-2,-l,-l,0 


0) 


2a9 4- ft + 2ft + ft + ft ft = 


= 




(-1,-1,-2,-2,-2,-1,0 


0) 

) 


4a9 - /3i + 2ft + /33 - ft + 2/?5 + ft + 2/?7 = 


= 


(_l,_l,_2,-3,-2,-l,0 


0) 


4Qg - /?i -h 2ft -h ft 3ft 2ft -1- /?6 


-f 2ft 


= 


(-1,-1,-2,-3,-2,-1,-1 


-1) 


4a9 + /3i + 2ft 3ft 3ft 2ft ft 


-2ft 


= 


(-1,-1,-2,-2,-2,-2,-1 


-1) 


4q9 -I- ^1 -I- 2ft -1- 3ft - ft 2ft ft 


■f 2ft 


= 


(_l^_l^_2,-3,-2,-2,-l 


-1) 


4a9 + ,81 + 2/32 3ft 3ft + 2ft, -f ft 


4- 2ft 


= 


(-1,-2,-2,-3,-3,-2,-1 


-1) 


4a9 + /3i -h 2ft 3ft - 3ft - 2ft - ft 


+ 2ft 


= 


(-2,-2,-3,-4,-3,-2,-1 




4a9 /3i 2ft - ft ft 2^5 + 3,flfi 


-f 2ft 


= 


(-2,-2,-3,-4,-3,-2,-1 


-1) 


4q;9 + /9i H- 2ft 4- 3ft -|- ft + 2ft 3ft 


+ 2ft 


= 


(-1,-2,-2,-4,-3,-3,-2 


-1) 


4a9 + 3;3i -f 2ft -f ft -f ft - 2ft - ft 


f 2ft 


= 


(-2,-2,-3,-4,-3,-2,-2 


-1) 


4Qg + 3/3i + 2ft + ft + ft -f 2/35 + 3ft 


-2ft 


= 


(-2,-2,-3,-4,-3,-3,-2 


-1) 


4q9 + 3/3i -f 2ft 4- ft ft 2/35 + 3ft 


+ 2ft 


= 


(-2,-2,-3,-4,-4,-3,-2 


-1) 


2a9 + /3i 2ft ft - ft ft = 


= 




(-2,-2,-3,-5,-4,-3,-2 


-1) 


2a9 /3i + 2ft -h ft + ft 4- ft = 


= 




(-2,-3,-3,-5,-4,-3,-2 


-1) 


4q;9 4- 3/3i 4- 2ft 4- ft - ft - 2ft ft 


f 2ft 


= 


(-2,-2,-4,-5,-4,-3,-2 


-1) 


2(19 + „5i + 2/32 + ft + ft + 2A^ + ft = 




(-2,-3,-4,-5,-4,-3,-2 


-1) 


4q9 4- 3/3i 4- 2ft -h ft - ft + 2ft + /36 


4- 2ft 


= 


(-2,-3,-4,-6,-4,-3,-2 


-1) 


4a9 + S^Si 2ft ft + 3ft + 2ft + ft 


+ 2ft 


= 


(-2,-3,-4,-6,-5,-3,-2 


-1) 


2a9 -h /3i -h 2ft ft ft - ft = 


= 




(-2,-3,-4,-6,-5,-4,-2 


-1) 


2a9 -h /3i -h 2ft ft ft ft = 


= 
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Case SAi 






o 1 


Pi = 


(0, 1, 0, 1, 1, 1, 


1, 


n^ 
0) 


o 2 


132 = 


= (1, 1, 2, 3, 3, 2, 


1 

1, 


n^ 


o 3 


= 


= (1, 1, 2, 2, 1, 1, 


1 

1, 


n^ 
Uj 


o 4 


Pi = 


= (0, 0, 0, 1, 0, 0, 


0, 


0) 


o 5 


= 


= (1, 0, 0, 0, 0, 0, 


0, 


n^ 
0) 


o 6 


Pe = 


= (1, 2, 2, 3, 2, 1, 


0, 


n^ 
0) 


o T 


P7 = 


= (0, 0, 0, 0, 0, 1, 


0, 


n^ 
0) 


o 8 


Ps = 


= (2, 3, 4, 6, 5, 4, 


Q 

3, 


9^ 


Ai-lowest 




Ijinear combination 










9r>,,^ R, M R^. ^ — 







f_l _l _l _2 -1 -1 0) 




9r»„ -1- fl. -1- flr -1- fl^. -1- fl-T — 

ZQg P4 P5 1 P6 P7 — 







f_l _i _2 -2 -1 -1 01 

\ -1-5 -1-5 ^5 -1-5 -i-,W5V^y 




9rvf^ R-, -\- R-^ -L R^ ~L Rt — 
Z09 — pi -f- p3 -f- p6 -t- p7 — 







f_l _1 _1 _2 -1 -1 -1 01 

\ -1-5 -I 5 -1-5 -^-i -1-5 -I 5 




9rv^ J- R. J- R^ -I- R, J- R^ — 
ZQg + Pl + P3 + P4 + Ps — 







f-l _i _2 -2 -2 -1 01 




zag — Pl P2 P4 P6 — 







-2-3,-2,-1,0,01 

\ -1-5 -1-5 iJ5 f-Tj 




9rvr, fl, -1- fl„ -1- fl, -i_ fl,. — 
ZQg — Pl + P2 + P4 + P6 — 







f_l _l _l _2 -2 -2 -1 01 

\ -1-5 -1-5 -1-5 ^5 ^5 




9rvn -1-/^1 -i- Rr, -i- R^ -I- R^ — 
ZCK9 -f- Pl -f- P2 P5 P7 — 







f_l _1 _2 -3 -2 -2 -1 01 

\ -1-5 -1-5 ^5 '-»5 ^5 ^^5 -i-5»-'y 




zag -f- P2 P3 P4 P7 — 







(-1,-1,-2,-3,-2,-1,-1,-1) 




ZQg — Pl + P4 ^ P7 + P8 — 







(-1,-1,-2,-2,-2,-2,-1,-1) 




9r>,,^ R^ J- R~^ J- R^ — 
zag — pi — P4 -f- P7 -f- P8 — 







(-1,-2,-2,-3,-2,-2,-1,0) 




9/^^ -I. R. M R^ M R^. R^ — 
ZQg + Pl + P3 + P6 + P7 — 







(-1,-1,-2,-3,-2,-2,-1,-1) 




9rvn R, 4^ R, 4- R^ 4- Ra — 
zag — pi P4 P7 P8 — 







(_l,_l,_2,-3,-3,-2,-l,-l) 




9/Ty^ R. 4- R^ R., 4- R,^ — 
zag — Pl + P2 ~ P3 + P8 — 







(-1,-2,-2,-4,-3,-2,-1,0) 




9rvn 4- R-, 4- Rr, 4- R^ 4- R^. — 
ZCKg -f- pi -f- P2 P4 P6 — 







(-1,-2,-2,-3,-3,-2,-1,-1) 




9or, flo Ra 4- Rc, 4- Ra — 
ZQg — P3 — P4 + P6 P8 — 







(-2,-2,-3,-4,-3,-2,-1,0) 




Orvr. 4- R^ 4- R., 4- R^ 4- R.. — 
ZQg + P2 + P3 + P5 + P6 — 







(-2,-2,-3,-4,-3,-2,-1,-1) 




9n.o — 4- /^e 4- flc 4- flo — 
zctg pi -f- p5 "T p6 -f- p8 — 







(-1,-2,-2,-4,-3,-3,-2,-1) 




2Qg + /3l + /34 + ^7 + ^8 = 







(-2,-2,-3,-4,-3,-2,-2,-1) 




2Qg + /?3 + ft - ;fl7 + ;fl8 = 







(-2,-2,-3,-4,-3,-3,-2,-1) 




2ag + /?3 + /?5 + /37 + ft = 







(-2,-2,-3,-4,-4,-3,-2,-1) 




2Qg + /32 - ft + ft + ft = 







(-2,-2,-3,-5,-4,-3,-2,-1) 




2Qg + ;92 + ft + ft + ft = 







(-2,-3,-3,-5,-4,-3,-2,-1) 




2Qg + /?1 + /?5 + ft + ft = 







(-2,-2,-4,-5,-4,-3,-2,-1) 




2Qg - /3l + ft + ft + ft = 







(-2,-3,-4,-5,-4,-3,-2,-1) 




2Qg + ft - ft + ^6 + ;fl8 = 







(-2,-3,-4,-6,-4,-3,-2,-1) 




2Qg + ft + ft + ft + ft = 







(-2,-3,-4,-6,-5,-3,-2,-1) 




2Qg + ft + ft - /37 + ft = 







(-2,-3,-4,-6,-5,-4,-2,-1) 




2Qg + /?2 + ft + /37 + ft = 







(-2,-3,-4,-6,-5,-4,-3,-1) 




2Qg + /?1 + ft + ft + ft = 







Let Ag-i be the root systems 


generated by all aj except j 


= i. 


Suppose that for some 


index i, 


the coefficient in front of ai is 


odd. 



Then, choosing a height function for which is Ag_i-lowest, we see 
that we land in one of the already treated cases (that of Ag-i = ^8, 
2A4, A2 + Eq 4j42 or As -f Ai -I- A2 with ai = 3) and we have a linear 
combination with a coefficient ±1 for some Qfc. On the other hand, 
the assumption that all coefficients in front of ai are even for all i 
immediately contradicts the minimality of the starting relation (|18|) . 

• A is of type F4,. All possibilities for Ai 7^ A come in seven conjugacy 
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classes ( |Dyn72| Table 10]). 





Case A2 + A2 






Bi = (1, 0, 0, 0) 




I. 


Bo = (0, 1, 0, 0) 




3 


B-i = (0, 0, 0, 1) 






04 = (2, 4, 3, 1) 




All possible as 


Linear combination 




(0,0,1,0) 


3^5 + + 4/32 + ft - ft 


= 


(0,1,1,0) 


3a5 + 2/3i + ft + ft - ft = 


= 


(0,2,1,0) 


3a5 + 2/3i - 2ft + ft - ft 


= 


(0,0,-1,0) 


3a5 - 2/3i - 4ft - ft -f ft 


= 


(0,0,1,1) 


3q5 + 2^1 + 4/32 - 2ft - ft 


= 


(1,1,1,0) 


3a5 - /3i + ft + ft - ft = 





(0,-1,-1,0) 


3a5 - 2/3i - ft - ft -f ft = 


= 


(0,1,1,1) 


3a5 + 2/3i + ft - 2ft - ft 


= 


(2,2,1,0) 


3a5 - 4/3i - 2ft + ft - ft 


= 


(0,2,1,1) 


3a5 2/3i - 2/32 - 2ft - ft 


= 


(0,-2,-1,0) 


3^5 - 2/3i + 2ft - ft + ft 


= 


(0,0,-1,-1) 


3q5 - 2/3i - 4/32 2ft + ft 


= 


(-1,-1,-1,0) 


3a5 /3i - ft - ft + ft = 





(1,2,1,0) 


3a5 - /3i - 2ft + ft - ft = 


= 


(1,1,1,1) 


3a5 - /3i + ft - 2ft - ft = 


= 


(0,-1,-1,-1) 


3^5 - 2/3i - ft 2ft + ft 


= 


(2,2,1,1) 


3a5 - 4/3i - 2/32 - 2ft - ft 


= 


(0,2,2,1) 


3a5 4/3i + 2ft - ft - 2ft 


= 


(-2,-2,-1,0) 


3^5 + 4/3i + 2ft - ft + ft 


= 


(0,-2,-1,-1) 


3q5 - 2/3i + 2ft 2ft -1- ft 


= 


(-1,-2,-1,0) 


3a5 + /3i + 2/32 - ft -f ft - 


= 


(-1,-1,-1,-1) 


3q5 + /3i - /32 + 2/33 -f ft - 


= 


(1,2,1,1) 


3^5 - /3i - 2ft - 2ft - ft 


= 


(2,2,2,1) 


3q5 - 2/3i + 2ft - ft - 2ft 


= 


(-2,-2,-1,-1) 


3q5 4/3i + 2/32 -f- 2ft -f ft 


= 


(0,-2,-2,-1) 


3a5 - 4/3i - 2/32 + Pi + 2ft 


= 


(-1,-2,-1,-1) 


Sas + ^1 -f 2ft + 2ft + ft 


= 


(1,2,2,1) 


3a5 2ft - ft - 2ft 


= 


(2,4,2,1) 


3q5 - 2/3i - 4ft - ft - 2ft 


= 


(-2,-2,-2,-1) 


3q5 4- 2/3i - 2ft ft -f 2ft 


= 


(-1,-2,-2,-1) 


Sas - /3i - 2ft + /33 2ft 


= 


(1,3,2,1) 


Sas /3i - /32 - /33 - 2ft = 


= 


(-2,-4,-2,-1) 


3a5 + 2/3i -1- 4/32 /33 + 2ft 


= 


(-1,-3,-2,-1) 


3a5 - /3i /32 -1- ft + 2ft = 


= 


(2,3,2,1) 


3a5 - 2/3i - ft - ft - 2ft 


= 


(-2,-3,-2,-1) 


Sag 2^1 + ft + /33 2ft 


= 


Number of relations: 36 





By direct observation we see that the mini-lemma from case E& 
applies to ft and ft for each of 36 possible values for as and we get 
that there exists some ai with coefficient = ±1. 
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Case Bi 




1 


/3i = (0, 1, 0, 0) 




2 


/32 = (0, 0, 1, 0) 




3 


/33 = (0, 0, 0, 1) 




, 4 


/34 = (2, 2, 1, 0) 


Possibilities for lowest 


as Linear combination 


(-1,-3,-2,-1) 


2a6 + 4j3i + 3/32 + 2/33 +134 =0 







Case Ai + A3 








/Ji = (l, 0, 0, 0) 








/32 = (0, 0, 1, 0) 




li 




/^3 = (0, 0, 0, 1) 








/34 = (2, 4, 2, 1) 




A 1 -lowest 


Q;b 


Linear combination 




(-1,-1,-1,- 


1) 


Aas +2l3i+ 202 +3133 +134 


= 


(-2,-2,-2,- 


1) 


2q.-, + 2/3i + 2/32 /33 /34 


= 


(-1,-3,-2,- 


1) 


4a5 - 2/3i 2/32 /33 3/34 


= 


(-2,-3,-2,- 


1) 


4q5 -h 2^1 2^2 +p3+ 3/34 


= 



Case C3+ Ai 

/3i = (l, 0, 0, 0) 

P2 = (0, 1, 

/33 = (0, 0, 

^4 = (2, 4, 



0, 0) 
0) 
2) 



A 1 -lowest as 

(-1,-1,-1,-1) 

(-2,-2,-1,-1) 
(-2,-2,-2,-1) 
(-2,-4,-2,-1) 
(-1,-3,-2,-1) 
(-2,-4,-3,-1) 
(-2,-3,-2,-1) 



Linear combination 
2a5 -2132 -133 +134 =0 

2a5 + 2/3i -133 + 134 = 
2a5 2/3i + /33 + /34 = 
2as + 201 +4:02+03+04 = 

2Q5 + 202 +03 +04 = 
2Q5 + 201 + 402 +303 + 04 =0 

2qb + 201 -h 2^2 -h /33 -h /34 = 





Case D4 






01 = {0, 0, 1, 0) 






02 = {0, 0, 0, 1) 






03 = {0, 2, 1, 0) 






04 = {2, 2, 1, 0) 




Ai-lowest as 


Linear combination 




(-1,-1,-1,-1) 


2a6 -h /3i -1- 2/32 -1- /34 = 




(-1,-2,-2,-1) 


2a5 + 201 + 202 +03 + 04 = 





(-1,-3,-2,-1) 


2a5 + 01 + 202 + 203 + 04 = 





(-2,-3,-2,-1) 


2aB +01 + 202 + 03 + 204 = 
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_ 1 


pi — (1, u, U, Uj 


I 


P2 — (U, 2., 1, UJ 


o 3 


P3 - (U, U, 1, UJ 


o -4 


_ ('9 A Q 9"! 

P4 — l,^, o, ZJ 


Ai-lowest Q5 


Linear combination 


(-1,-1,-1,0) 


2q5 + 2/3i /32 ft = 


(-1,-1,-1,-1) 


2^5 - ^2 -1- ft = 


(-2,-2,-1,-1) 


2q5 + 2/3i - ft + ft = 


(-2,-2,-2,-1) 


2a5 2/3i ft -f ft = 


(-1,-2,-2,-1) 


2a5 ft + ft = 


(-2,-4,-2,-1) 


2a5 + 2/3i -f 2ft - ft + ft = 


(-1,-3,-2,-1) 


2az, + ft + ft = 


(-2,-4,-3,-1) 


2^5 + 2/?i + 2ft + ft + ft = 


(-2,-3,-2,-1) 


2q5 + 2/3i + ft + ft = 





Case 4yli 




o 1 


ft = (2, 2, 1, 0) 




o 2 


ft = (0, 2, 1, 0) 




o 3 


ft = (0, 0, 1, 0) 




o 4 


ft = (2, 4, 3, 2) 




Ai-lowest as 


Linear combination 




(0,-1,-1,0) 


2^5 -1- ft -f ft = 




(-1,-1,-1,0) 


2a5 -1- ft -f ft = 




(-1,-2,-1,0) 


2a5 -f ft + ft = 




(-1,-1,-1,-1) 


2^5 - ft + ft = 




(-2,-2,-1,-1) 


2a5 -1- /?i - ft - ft + ft = 





(-2,-2,-2,-1) 


2a5 -f /?i - ft + ft + ft = 





(-1,-2,-2,-1) 


2^5 + ft + ft = 




(-2,-4,-2,-1) 


2^5 /?i + ft - ft + ft = 





(-1,-3,-2,-1) 


2^5 + ft -1- ft = 




(-2,-4,-3,-1) 


2a5 -1- /3i ft ft ft = 





(-2,-3,-2,-1) 


2^5 +Pi+Pi = 





If the coefficient as of as is odd then Ai = A2+A2 or F4, and we get 
the desired relation with coefficient at — ±1 for i 7^ 5 in the former 
case, and as = 1 in the latter. On the other hand, the assumption 
that all coefficients in front of a; are even for all i immediately 
contradicts the minimality of the starting relation (|18p . 

• A is of type G2. If the three roots ai are of equal length and positive 
then one of them is the sum of the other two. If two of the at are 
long, say ai and a2, and the third, as, is short, then by direct check, 
after eventually switching ai and a2, we get that 3a3 = ±2ai ± 02- 
If two of the ai are short then the third root is their integral linear 
combination. The general case for arbitrary positivity of the roots 
ai's follows directly. 
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6 On the actual implementation 



6.1 Modifications to Section 12.31 used in the im- 
plementation 

In this section will stand for the root system of a classical Lie algebra 
and will be extended to a set A. In addition to the data given for the 
algorithm in Section [2.31 we will give a total order on A for which a\ >- 
• • • ^ ajv. In the context of computing the vector partition function, the 
total orders >- were defined in (BBCV OG . 

We are using the notation from Section [1.1. II For define A to be 
the root system of A„, define -< to be the order induced by rii — r]i+\ >- 
Tji+i — rji+2 and by a + /? ^ ct, P- Equivalently, define >- to be the order 
naturally induced from the lexicographic order on weights rji of gl(n), 
projected down to sl(n). For systems D„, respectively B„, define A to be 
the root system of B„, respectively D„, U{2r;i} and define -< in the same 
way as for A„. For system Cn define A to be the root system of C„ and 
define the order in the same way as for A„. 

6.1.1 An algorithm for the partial fraction decomposition 
of the generating function for the vector partition function 
of the classical root systems 

Step 1 (Reduction step) . We are given a fraction p{x) — ^ . . . -^-^ — i-^ , 
ii < • ■ • < ifc . 

Step 1.1 Find the smallest index a, and the smallest index b depending 
on a, and the smaller of the numbers n = 1,2, for which — 
nai^ G A and a;^ — nai^ -< ai^. If no such couple of indices 
{ia,ib) exists label the fraction as "reduced". 

Step 1.2 For the couple of indices {ia,ib) perform the reduction formula 
Lemma I2.1[ e'). Label any newly created fractions as "non- 
reduced" . 

Step 2 (Termination step). Find the first fraction labeled "non-reduced" in 
the memory and go to Step 1. If all fractions are labeled as reduced 
proceed to Step 3. 

Step 3 (Getting rid of the redundant short roots in the denominators). 
Transform the output of Step 2 to a partial fraction decomposi- 
tion (see Definition I2.3|l . One needs only to apply the formula 
1^7^ = i-x'^c where it is needed (see Lemma |6. II below) . 

We leave the following lemma without prooQ. 

Lemma 6.1 Let A"*" C A 6e the set of positive roots of a classical root 
systems of rank n. Steps 1 and 2 of the algorithm in Section \6. 1 . 1\ vroduce 
a decomposition of the generating function of the partition function of any 

*The "vector partition function" program is designed to display an error message and crash 
if it doesn't reach a partial fraction decomposition. 
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subset I C A"*" of the form 

n 2 _ — X/ Ph,...,l„,mi,...,m„,cxi^,...,ai„{x) 



- ■ X' 



1 



(1 - x«i )'i (1 - ■■■(!- x°''n )'»(!- a;2«in ' 

where 71 is exactly the VGfik of the voot systetu and Pii,...,in,mi ,...,mn jCKii 
are polynomials and Xi ... .x„ ■ 



6.2 Short outline of program functionality 

The "vector partition function" program is being written in C++. It 
uses template classes for polynomials over arbitrary commutative rings 
with identity, a realization of the ring of quasi-numbers (multiplicative 
relations realized, additive not), classes for operations with large rational 
numbers (size limit dictated only by system resources), and hash-table 
array management template classes. Additionally, there are classes for 
polyhedral cone divisions of Euclidean space, classes for computing orbits 
of Weyl groups and Kazhdan-Lusztig polynomials and coefficients. The 
entire code is written from scratch by the author; no external packages 
are used. The standard C++ library for stream/string manipulations is 
used only for the final printouts of the algorithm. 

Before and after each partial fraction decomposition runs, the "vector 
partition function" program is set to substitute rational numbers in the 
expressions it transforms. Comparing the two numbers has been our main 
tool for catching programming mistakes. Each step of the program has 
successfully passed such "checksum" tests. 
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A Appendix. Complex exponents and 
quasi- numbers 

Let X €Z. Then 

= l^:^-"'"' = T n (e"*-e"*) 

0<fc<i3 
k^l 
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B Appendix. Partial fraction decompo- 
sition tables for A2, A3, A4, B2, B^, C2, C3 
and G2 

Ill 1 1 1 

(1 — Xl) (1 — X2) (1 — XlX^) (1 — X1X2) (1 — xfx^) (1 — X1X2) 



+{xiX2 + X1X2 * + X1X2 ^ + 4x^X2 ^ + '^x\x2 + ^x\x2 ^ + 6xja;2 + 6x10:2 ^ + &x\x2 ^ 

+4x^X2^ + 4xtx2^ + 4x\x2^ + x\x2'^ + x\x2' + x'{X2^) 

1 1 

{l-xlY (l-xixa) 

~t~(l3X2 ~t~ 15^2 ~t~ Xj X2 ~t~ X2 ~t~ 2Xj^ 0^2 ~t~ 2X]^ X2 ~t~ 5X;^ ^2 ~1~ 6X]^ X2 ~t~ 6X;^ X2 

~t~9iCi 2^2 ~l~ llcC^ X2 ~t~ Xj iC2 ~l~ X^ X2 ~\~ X^ X2 ~t~ ScC]^ X2 ~\~ X2 ~\~ 9X]^ X2 

, -2 -13 , c -1 -7 , o -1 -6 , -2 -7 , -2 -8 , -1 -4 , -1 -5 , -3 , r -4 

+11x1 a;2 + 0x1 X2 + X2 + Xl X2 + Xl X2 + Xi X2 + .x'l X2 + X2 + 0x2 
+5x2 ^ + 2xix^^ + 6x1X2 ^ + 11x1X2 + X1X2""'" + x'i + 5x^X2^ + 9xix2^ + X1X2 + 5x'i + IIX1X2 
+12x1X2"* + 10x1X2 + llx^-^x^" + 14x^^° + 4x1^X2'^'' + 10x3:^x2 " + 16x3:^X2 + 14X]"^x^^^ 

+ 10X2^'^ + 8X2 ^ + 10X1X2 + 14x^X2^ + 14X1X2* + 16x'iX2 ^ + 4xiX2 + lOXi + 14x\x2^ + 3xiX2 * 

-x\x2^ + 16x^" + llxix^^° + 6xr^X2 + 10x^^X2 " + 14x2 + lla;^^* + lOxix^^ + 12xiX2 ® 

+10x1X2"'' + 16xjx2 ^ + 11x^X2'' + 14XiX2 ^ + 6xfx2 + lOxf + llx\x2^ — 2x?X2 * — 3x?x^® 
+9x1X2""'"' + 5x1X2 ■'"'' + 4xi^X2^^ + 5x2^"^ + 6x1X2 + 5x1X2""""^ + 5x^X2'^ + 8x?X2 ^ + 5x?X2 ^ 

+9X1X2 ^ + 5X1X2 + dxtX2^ + 4XiX2 + 5xf + 5X1X2 "^ — X1X2 * + 2x\x2^^ + X?X2 ""^ 

+X2 + X1X2 + X1X2 + X1X2 + X1X2 + 2X1X2 + X1X2 + 2X1X2 + X1X2 
+X?X^^ + X1X2 + xl + xjx2 ^) 

1 1 

(1-xir (i-xixi) 
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(1 - (1 - X1X2) 

+(— 24.T1X2 * — 15x1X2 ^ — 8x1x2 ^ — 3x1X2 ^ — 8x1X2 * — llxix^^ — 13xixJ^ — 2x^^ — 3a;2 

e -7 r -8 -1 -7 -1 -8 -1-9 2 r, 2 -1 3 3 -, n 3 -1 

—0X2 — 0X2 — Xi X2 — Xi X2 — x'l X2 — x'l — 2x1X2 — X1X2 — 6x1 — 10x1X2 
—x\xi — 5xtx2 — 15xi — X1X2 ^ — 23xf X2 ^ — X2 — 36x^X2'' — 30xfx2 ^ — 24xf xj^ 
— 22xfx2 ^ ~ 14a;ia;^^ — 10x1X2 * — 4x2 ^ ~ 20xix^^ — 4xixi — 10xfx2 — 20xf — 32xiX2 ^ 
— 29x1X2 "* — 30xiX^^ — 26x1X2"' — 18x1X2 ^ — 16x1X2^ — lOxiXj *^ — 6x1X2 ^ — 20x1X2 ^ 
—6x1X2 — lOxf X2 — 15xi — 23xiX^^ — 12xiX^* — 15xiX^^ — IAX1X2'' — 7xfx2^ 

n3-7 r3-8 ,2-9 in 6-1 j72 r7 c7 06 -2 r,6-4 

— 9X1X2 — 0X1X2 — 4X1X2 — 10X1X2 — 4X1X2 — 0X1X2 — 6X1 — 8X1X2 — 2X1X2 

06-3 06-6 6-6 r) 4 -7 4 ,-8 ,3-9 r, 7 -1 8 2 
— 0x1x2 — 0x1x2 — x\x2 — zx1x2 — x\x2 — x1x2 — ^x\x2 — x1x2 

8 8 7 -2n 

— X1X2 — Xi — X1X2 ) 

1 1 



(1-X2)5 (1-xfxi) 



, I -4 -13 -4 -14 -4 -IS r> -3 -12 a -3 -13 r- -3 -14 „ -2 -11 

+(— Xi X2 — Xi X2 — Xi X2 — 2xi X2 — 6x1 X2 — 0x1 X2 — 7xi X2 

— 10X1^^X2""'^^ — 15X]^^X2 "'^ — X~[^X2^^ — 2X1^X2^'^ — 6X1^X2^ — 10xj~"'^X2 "'^ — 19x^^X2^^ 

— X 3^2 — — "^2 — ^"^2 — 153^2 — 43^-]^ ' X'2 — 103^]^ •^'2 — 20^^]^ 3^2 203^]^ 3^2 

^2^" — 26x2^^ " 4xix^^ — lOxix^* — 20xix^^ — 6x^^xJ"'^^ — lOx^^x^"*^* — 20x^"'^^ — 15x2 
;iX2 ^° — 19xiX2^^^ — 6x1X2 ^ — 10x1X2 * — 15x1X2'^ — 4x1^X2^'' — 5x2 ^* ~ 10x1X2 

CC^ 10cC]^iC2 TX\Xi2 4cC]^3/2 ^X'^X'^ 6iC]^CC2 2^2 X\X^ ^Xi^X'2 

-13 r, 3 -10 3 -11 4 -7 4 -8 4 -9n 
■,2 — 2X1X2 — X1X2 — X1X2 — X1X2 — XiXo ) 



-20X2 ^° - 
— 20X1X2 ^° 

—6X1X2^^— 1^ 1^ - 

2 -13 r, 3 -10 3 -11 4 -7 4 -8 4 - 
— X1X2 — 2X1X2 — X1X2 — X1X2 — X1X2 — X1X2 

1 1 

(1-X?)5 (1-Xlxi) 



+ (X2 ® + 4xiX2 ® + 6X?X2 ® + 4x?X2 ^ + xtx2 ^) 



(1-X?)5 (I-X2) 



, / 3 , 4 , 5 2 I 2 -2 I 3 -1 I 4 , -4 , 2 -3 , 3-2 

+ (Xi + X1X2 + X1X2 + X1X2 + X1X2 + Xi + X1X2 + X1X2 + X1X2 



(1-Xl)5 (1-xfxi) 

I /■ -4 , 2 -3 -4 -3 2 -2 2 -1 3 4 n 1 1 

+ (XlX2 + X1X2 — Xj — X1X2 — X1X2 — X1X2 — Xi — X1X2) 



'(1-Xl)8 (1-Xlxi) 



1 

= -X2- 



, ^ X« (1 - X1X2) (1 - X2)2 

1 1 



3 

-3:2X3- 



(1-Xl) (1-X2)2 

1 11 



^^^^ - ■ X" (1 — X1X2X3) (1 — X2X3) (1 — Xs)^ 

2 1 1 1 

-X2X3 



+2x1 



(1 - X1X2) (1 - X2X3) (1 - X3)'' 
111 
(1 -xi) (1 -X2X3) (1 -xs)-* 



2 1 1 1 

+3:2X3- 



3 

-X2X3 



+ 



(1 — X1X2X3) (1 — X2) (1 — X3)* 
1 11 

(1 — X1X2X3) (1 — X2X3)^ (1 — X3)^ 
111 



(1-Xl) (1-X2)2 (1-X3)3 
111 



-X2 



+X2X3 
-2X3 



(I-X1X2) (1-X2)2 (1-X3)3 
1 1 1 



(I-X1X2) (I-X2) (I-X3)* 
1 1 1 



1 2 
+3:3 



(1-Xl) (I-X2) (1-X3)4 



(1-Xl) (1-X2X3)2 (1-X3)3 
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_|_2 3 6-'- ^ 

(1 — X1X2X3X4) (1 — X2X3X4) (1 — X3X4) (1 — X4y 

3 5 1 1 1 1 

— 2a;2a;3a;4 



(1 — a;ia;2a;3) (1 — X2X3X4) (1 — X3X4) (1 — 0:4)'' 



+4x2X3X4 



(1 — a;ia;2) (1 — a;2a;3a;4) (1 — a;3a;4) (1 — X4y 

—4x3X4- 



(1 — xi) (1 — a;2a;3a;4) (1 — X3X4) (1 — X4y 



+2X2X3xi 



3 6 

+a;2a;3a;4 



(1 - X1X2X3X4) (1 - a;2a;3) (1 - a;3a;4) (1 - X4)'' 



(1 — X1X2X3X4) (1 — a;2a;3a;4) (1 — a;3a;4)'' (1 — X4Y 



_4^^^2^5 1 1_ 1 1 

(1 — a;ia;2a;3a;4) (1 — X2) (1 — X3a;4) (1 — X4y 

-2x2x1x1 ^ ^— 1 1 

(1 — X1X2X3X4) (1 — a;2a;3a;4) (1 — 0:3) (1 — X4y 

+2X2X3X1- 



„3.6 1 111 



(1 — X\X2X3X4) (1 — a;2a;3a;4) (1 — a;3a;4)3 (1 — X4Y 



+2X2X%x\ 



(1 — X1X2X3X4) (1 — X2X3X4) (1 — X3X4y (1 — X4Y 

2 3 1 1 1 1 

(1 — xiX2) (1 — X2a;3X4) (1 — X3y (1 — 0:4)'^ 

2 2 1 1 1 1 

+X3X4- 



(1 - Xl) (1 - X2X3y (1 - Xs) (1 - X4Y 

— 2a;2a;3a;4 



(1 - X1X2X3) (1 - X2a;3) (1 - X3X4) (1 - X4y 

^^^^^2^4 1 1_ 1 1 

^ * (1 — X1X2X3) (1 — a:2) (1 — a;3a;4) (1 — X4y 

^ 2 4 1 1 1 1 

+a;2a;3X4 



(1 — X1X2X3) (1 — X2) (1 — a;3a;4)^ (1 — X4y 



3 3 
—X2X3X4 



(1 — X1X2X3) (1 — a;2a;3a;4) (1 — xsy (1 — X4y 



43 



(1 — X1X2X3) (1 — X2X3X4) (1 — X3) (1 — X4y 
^^^3^5 1 1 1 1 

(1 — X1X2X3) (1 — X2X3X4) (1 — X-iXi)^ (1 — Xi)^ 
o 2^5 1 1 1 1 

(1 - xi) (1 - X2X3X4) (1 - a;3a;4)'' (1 - 0:4)* 

_4^^^2^4 ^ 1 1^ 1 

(1 — X1X2) (1 — X2X3) (1 — a;3X4) (1 — X4y 

3 1 1 1 1 

+0x2X4- 



2 

-X3X4 



(1 - 0:10:2) (1 - 0,-2)^ (1 - X3X4) (1 - 0:4)'^ 



(1 - xi) (1 - a;2a;3)2 (1 - 0:3)2 (1 - a;4)s 



—Ax2x\x\ 



1 2 5 

+0:20:3X4 



(1 — 0:10:2) (1 — 0:20:30:4) (1 — 0:3) (1 — 0:4)'' 



(1 — 0:1X2) (1 — X2X3X4) (1 — 0:30:4)* (1 — 0:4)" 



+2X2X3x1 



(1 — 0:10:2) (1 — 0:20:30:4) (1 — 0:30:4)^ (1 — 0:4)^ 

^ 3^^^2^5 1 1 1 1 

(1 — 0:10:2) (1 — 0:20:30:4) (1 — 0:30:4)2 (1 — X4)® 

~\~^x'^x^ 

(1 — 0:1) (1 — 0:20:3) (1 — 0:30:4) (1 — 0:4)'' 

2 5 1 1 1 1 

-0:3X4- 



(1 — Xl) (1 — X2X3X4)2 (1 — X3X4) (1 — X4)'= 
2 1 1 11 

"^^^ (1 - Xl) (1 - X2X3)2 (1 - X3)3 (1 - X4)4 

~\~^x^x^ 

(1 — Xl) (1 — X2X3X4) (1 — xa) (1 — X4)'' 

^^2^5 1 1 1 1 

(1 — Xl) (1 — X2X3X4) (1 — X3X4)^ (1 — X4)^ 

-4x1x1^ ^ ^ ^— 

(1 — Xl) (1 — X2X3X4) (1 — X3X4)2 (1 — X4)® 



44 



^ " (1 — X1X2X3X4) (1 — X2X3) (1 — Xsy (1 — X4)^ 

3^.^^3^4 1 1 1 1 

^ (1 - XxX2X3X^ (1 - X2Xz) (1 - Xz) (1 - X^^ 

3 5 1 1 1 1 

-\-X2XzX4 



2 5 
-X2XZX4 



, 3 6 

+0:2X314 



(1 — X1X2X3X4) (1 — a;2a;3) (1 — xzX4Y (1 — X4Y 



(1 - x-ix-2X3a;4) (1 — X2) (1 - X3X4)'' (1 — X4)'' 



(1 — X1X2X3X4) (1 — a;2a;3a;4)^ (1 — a;3a;4) (1 — X4Y 



+Ax2X%x\- 



(1 — X1X2X3X4) (1 — X2) (1 — 0:3) (1 — X4y 



—2x2x'ix\ 



-3x2X3X4 



(1 — XlX2XiX4) (1 — X2) (1 — X3a;4)^ (1 — X4)^ 

1 111 

(1 — X1X2X3X4) (1 — X2) (1 — a;3a;4)^ (1 — X4Y 



"^'^^ (1 - Xl) (1 - X2) (1 - X3)2 (1 - X4)6 
2 1 1 1 

"^'"^^ (1 - X1X2) (1 - X2X3) (1 - X3)4 (1 - X4)4 
1111 



+3X2X4 



(1 - X1X2) (1 - X2)2 (1 - X3)2 (1 - X4)5 



<j 2 1 11 1 

^^'^''''''^ (I-X1X2) (1^^ (1 - X3)2 (1-X4)6 

—4X3X4 



+4X3X4 



(1 - Xl) (1 - X2X3) (1 - X3)3 (1 - X4)5 

1111 
(1 - Xl) (1 - X2) (1 - X3)3 (1 - X4)5 



^ 2 1 1 1 

+2x|- 



' (1 - Xl) (1 - X2X3) (1 - Xs)" (1 - X4Y 



45 



2 1 111 

"^"^'"^^ {I-X1X2X3) (1 - X2) (1 - XS^ (1 - X4r 

^ 3 3 1 1 1 1 

+X2X3X4 ^ :7-7T9 r-TT ZTTe 



3 



—X2X3 



(1 — X1X2X3) (1 — X2X3)^ (1 — X3X4,) (1 — X4)'^ 
1 111 



4 



(1 - X1X2X3) (1 - a;2a;3) (1 - X3Y (1 - X4) 

_^3_^ 1 1 1 1 

(1 — X1X2X3) (1 — X2X3) (1 — Xs)^ (1 — Xa)^ 

3 2 1 1 1 1 

(1 — X1X2X3) (1 — a;2a;3) (1 — xs)^ (1 — xa)'' 

^ o^^^3^3 1 1 1 1 

(1 - X1X2X3) (1 - a;2a;3) (1 - X3) (1 - a;4)'' 

2 1 111 

(1 — xia;2a;3) (1 — X2) (1 — xs)^ (1 — a;4)^ 

2 2 1 1 1 1 

+3x2x3x4 _ ^^^^^^^ _ _ _ 

-4x xlx^ ^ 111 

(1 — X1X2X3) (1 — X2) (1 — X3) (1 — X4)'' 

_^ ^3^2 1 1 1 1 

(1 — X1X2X3) (1 — X2X3)^ (1 — X3) (1 — X4)® 
^ 3 1 1 1 1 

^(1-Xl) (1-X2)2 (I-X3X4) (1-X4)6 

4 1 1 1 1 



-2X2X3X4 



+2X2X3X4 



2 2 
-8X2X3X4 



+4X2X3X4 



(1 - Xl) (1 - X2) (1 - X3X4)4 (1 - X4)4 
1111 



(1 


- X1X2) (1 


-X2) (1-X3)3 (1- 


X4)5 




1 


1 1 


1 


(1 


- X1X2) (1 


- X2X3) (1 - x-sY'' (1 


-X4)5 




1 


1 1 


1 


(1 


— X1X2) (1 


— X2X3) (1 — xs)^ (1 


- X4)'^ 




1 


1 1 


1 



(1 - X1X2) (1 - X2X3) (1 - X3) (1 - X4)'' 



46 



2 4 1 1 1 



6 



(1 — X\X2) (1 — X2Xz) (1 — XzXiY (1 — Xi)^ 
_g^^^2 1 1__ 1 1 

" (1 - X1X2) (1 - X2Y (1 - X3) (1 - X4) 
3 1 1 1 1 

(l-a;ia;2) (1 - 0:2)2 (1 - x^x^)^ (1 - x^Y 
^3 1 1 1 1 

(1 — X1X2) (1 — X2Y (1 — XzXiY (1 — XiY 

o 3 1 1 1 1 

—6Xi- 



- (1 - Xl) (1 - X2)2 (1 - X-yX^Y (1 - 0:4) 



5 



— a;2X3a;4 



1111 



(1 — X1X2) (1 — X2) (1 — a::3a^4)'* (1 — xY) 

O 4 1 1 1 1 

— 2X2X3X4^ 



4 



— 2X2X3x\ 



(1 — a;ia;2) (1 — X2) (1 — a;3a;4)^ (1 — xY)^ 
1111 

(1 — X]_X2) (1 — X2) (1 — xsxY'^ (1 — a;4)'' 



. 2 2 1 1 1 1 

+ ^33:4 _ _ ^^^^^ _ ^^^2 _ ^^^6 

,231 1 1 1 



—x\ 



(1 - Xl) (1 - a;2a;3) (1 - xs) (1 - xY)'' 
1111 



4 



(1 - xY (1 - 2^2)2 (1 - a:3a;4)3 (1 - xY 
, 2 1 1 1 

+XzX4 _ _ X2X3XY^ (1 - 0:3) (1 - xYl^ 
^2^5 1 1 1 1 

(1 — Xl) (1 — a;2a;3a;4)2 (1 — a;3a;4)* (1 — x^Y 
2 5 1 1 1 1 

-XzXi- 



' (1 — xY (1 - X2XzxY'^ (1 — XzxY'^ (1 - xY'' 
„ 2 1 1 1 1 

(1 - Xl) (1 - 0:2)2 (1 - 0:3) (1 - xaY 

4 1 1 1 1 

+ X3X4 _ _ _ ^^^^^3 _ ^^^5 



4 1 1 1 1 

+ X3X4 _ _ _ ^^^^^2 _ ^^^g 

^.^^.3^,5 1 1 1 1 

(1 — X1X2X3X4) (1 — X2X3X4)2 (1 — X3) (1 — X4)'^ 

3 6 1 1 1 1 

(1 — X1X2X3X4) (1 — X2X3X4)2 (1 — X3X4)3 (1 — X4)'* 

~|~tZ'2«i'3 "-C^ 

(1 — X1X2X3X4) (1 — X2X3X4)2 (1 — 13X4)2 (1 — X4)^ 

1111 



-3X4- 



(1 - Xl) (1 - X2)2 (1 - X3)2 (1 - X4)5 

1 1 1 1 

(1 - Xl) (1 - X2)2 (1 - X3)3 (1 - X4)4 
1111 



(1 - X1X2) (1 - X2)2 (1 - XzY (1 - X4)'' 
1111 



-I-X2X3 
-2X3 



(1 - X1X2) (1 - X2) (1 - X3)4 (1 - X4)4 
1111 



(1 - Xl) (1 - X2) (1 - X3)'' (1 - X4)4 

2 3 1 1 1 1 

-X3X4- 



(1 - Xl) (1 - X2X3)2 (1 - X3X4) (1 - X4)^ 



3 

-X2X3 



(1 — X1X2X3) (1 — X2X3)2 (1 — XzY (1 — X4)'' 



+X2xix4 ' 



(1 — X1X2X3) (1 — X2XzY (1 — (1 ~ 3:4)® 



47 



— X2 



(1 - XxX-i) (1 - 



+ 



1 1 



(1-Xl) (1-X2)3 



n 1 - a:° ^ +(-a;?xi - 2a;?a;2 - xixi) 



aSC2 



(l-a:?X2) (1-2:2)3 



2 1 
+a;2 



(1 - x\xi) (1 - a;2)3 

1 1 



+(-a;i + l) 



(1-Xi)2 (I-X2) 



2 



.2 X 1 1 

+(a;ia;2 + x\Xi - xi)- 



(l-xi) (1- 0:2)3 



nl , / 2 12 , 3 14 11 o 2 11 , o 3 13 o 10 o 2 10 , o 3 12 o 9 2 9 

ot ^ +{-a:2a;3 + 2^2X3 - - ^x^^x^, + 3x2X3 - 3x3 - 3x2X3 + 3x2X3 - 3x3 - X2X3 

, 3 11 8\ 

+X2X3 -X3) 

1 1 1 

(1-Xl)(l-X2xi)2 (1-Xi)6 

+ (X2X- + 3X2xf + 3X2XJ^ + X2X-) ^^-^^ (T^^^ (T^^ 

+ (-X2xJ° - 6X2X3 - 15X2X3 - 20X2X3 - 15X2X3 - 6X2X3 - X2X3)-- i --; 3^ --; 

(1 - X1.T2) (1 - X2X3) (1 - X5) 

+ (X3" + X2X3^ + 10X2X3^ + 10X2X3^ + 4X3 + 4X2X3" + 40X2X3^ + 50X2X3^ + X2X3 + 6X2X3 

+63xix3° + 105xix3° - 18x3 + 4x2x1 + 49xix| + 119x1x1 - 30x| + X2X3 + 19xlx| + nx\x% 

, n 7 8 , -I -I 7 9 , o 6 7 I n-7 5 7 I r, 2 11 r, 3 13 , r 7 12 , , , 2 10 
+2X2X3 + 11X2X3 + 3X2X3 + 27X2X3 + 2X2X3 - 2X2X3 + 5X2X3 + 11X2X3 

-llxixg^ + 17x^x" + 21xix| - 21xix" + 21x^xJ° + 17xix| - 17x1x3° + 5xix| 

1- 3 9 4 14 , 8 12 o 4 13 , o 8 11 o 4 12 , „ 8 10 4 11 
— 5X2X3 — X2X3 + X2X3 — 3X2X3 + 3X2X3 — 3X2X3 + 3X2X3 — X2X3 

-21x| - 7xt + 4xlxt - xl) 

1 1 1 

(1-Xi) (I-X2) (l-xi)7 

+(—4x3 ~ a;2x| — 6x2X3° — 3x3 — 4x2x| — 26x2X3 + 6x3 — 6x2X3 — 45x2x| + IIX3 — 4x2x| 
—39x2X3 + 6x1 — 3^2X3 — 17x2xf + X3 — 10x2X3 — 3x2X3 — X3 — 2x2x1 + 2x2x" 

^6 10 in28,in3 10 ,,6 9 io27|,o39 1068 i^26 
—4X2X3 — 10X2X3 + 10X2X3 — 14X2X3 — 18X2X3 + 18X2X3 — 18X2X3 — 14X2X3 

+14xix| — 4xix| + 4xix3 — 2x2x| + x|x3^ — X2X3° + 3xtx" — 3X2x| 

+3X2X3° - 3X2X3 + X2xl - X2X3) 
111 

(i-xi)(i-x2)2 (i-xir 



48 



+(-xlxl^ - -ixlxf - Qxlxf - -ixlxf - 2xlxf - <dxlxl^ - IQxlx^^ ~ Mxax" 

— X%X^ — QX2x]^ — lbx%x]^ — 20X2X3" — x\x% — &X%X% — 15X23^3 — 3:2X3 

-6x%x% - x\x\) 

1 1 1 

(1-Xl) (1-X2X3)2 {1-xlY 

+ (4X2X3* + 14X2X" + 20X2X3^ + 15X2X" + 6X2X3° + X2X3 + 6X2X3"* + 16X2X3^ 

+ 15X2X3^ + 6X2X" + X2X3° + 4X2X3* + 9X2X3^ + 6X2X3^ + X2X" + X2°X3* 

, r, 10 13 , 10 12n 
+2X2 X3 +X2 X3 ) 

1 1 1 

(1 - Xix|x|) (1 - X2X3)^ (1 - X%Y 

+ (-X2X3° + X2X3 - 3X2x1 + 3X2X3 - 3X2X| + 8X2X3 - X2X3 + X2X3) 
1 1 1 

(1 - xixixi) (1 - X2Y (1 - 

+(2xix3^ - 5x1x3* - 10x|x3* - 10x2X3* + llxixg* - ITxIxg^ - 40xix" - SOxjx" 

I ni 5 13 01 9 12 8 12 -, 7 12 . -, - 5 12 , , 9 11 8 11 iin 7 l: 

+21X2X3 — 21X2X3 — 63X2X3 — 100X2X3 + 17X2X3 — 11X2X3 — 49X2X3 — 119X2X3 

, r 5 11 n 9 10 in 8 10 7 10 , 6 16 10 14 , o 6 15 o 10 13 

+0X2X3 — 2X2X3 — 19X2X3 — 77X2X3 + X2X3 — X2 X3 + 0X2X3 — 3X2 Xz 

+3X2X3* - 3X2"X3^ + X2X3* - X2"x" - 27X2X| - 4X2X3 - 3X2X3) 

1 11 
(1-Xlxixi) (1-X2)(l-Xi)^ 

+ (-2x2x" + 4x2xf + 6x2xJ^ - IQxlxf + 14x2x" + 26x2x" - ISxax" + 18x3X3° 
+45x^X3° - 14x^X3° + lOxixl + 39x^x1 - 4x|x| + 2x^x1 + 17x^x| - xlxg* 

+X2X3^ - 3X2X3^ + 3X2X3* - 3X2X3^ + 3X2X3° - X2X3* + X2xl + 3X2X3) 
1 11 

(1-xixixi) (1 - X2)2 (1 - a^ir 



+ (2X2X3* — 3X2X3° + 9X2X3° — 11X2X3 + 15X2x| — 15X2X3 + llX2x| — 9X2X3 

,o67 r>76, 6 12 8 10, 06 11 o89|o6 10 qBS 
+3X2X3 — 2X2X3 + X2X3 — X2X3 + 3X2X3 — 3X2X3 + 3X2X3 — 3X2X3 

+xix| - xlxs) 

1 1 1 

(1 - xixixl) (1 - X2)a (1 - xg)5 
+(-xix3° - 3x^X3'^ - 3xixJ* - xix"^ 



(1-Xlxixi) (1-X2xi)2 (1-Xi)6 

+(-10x2X3^ + 7x3° + x" - 50x2X3^ + 21x3 - 105x2X3* + 35x3 - 119x2X3° + 35x3 - 77x2xi 
+21x1 - 27x2x| + 7xi - 4x2X3 + X3 - 10x2X3* - 40x2X3^ - 63x2X3* - 49x2X3° 

iri69 068 r7 13 -,»7 12 r>i7 11 ii7 10 r>79 8 13 
— 19X2X3 — 3X2X3 — 0X2X3 — 17X2X3 — 21X2X3 — 11X2X3 — 2X2X3 — X2X3 

—3x1x3^ — 3x1x3* — 3:2X3°) 



(1-xi) (I-X2X3) (l-x§)7 

+ (X2X| - xlxl + 3X2X3 - 3X2X1 + 3X2X| - 3X2X3 + X2xl - X2X3) — — ^ 

(1-xi) (1-X2)4 (l-a;g) 
+(-2xixJ* + xixj° + 16x2x| - llxlxf + 4xix| + 7x2x| - 21xix| + 6xix| + 18x2x1 

-17xix| + 4x1x3 + 23X2x1 - 5x1x3 + xlxl + X2xl + 6X2X3 + X2xl) 

111 

(1-X1X2) (1-X2) (i-xi)7 

+(2x1x1 — x%x% — 6x2x1 + 10xlx| — 4x1x3 — 6x2x1 + 18x1x3 — 6x1x1 — 13x2x1 
+14x1x1 — 4xlx| — 13x2x1 + 4x1x1 — xlxl — X2X3 — X2xl) 

1 1 1 

(I-X1X2) (1-X2)2 (1-X^)« 



49 



+ {—2x2x1 + x\x% + X2X% — 'dx%x% + 4:X2xl + 5x2X3 — 15x2X3 + 6x^X3 + 9x2x1 — 11x2X3 
+4x1x3 + 7x2X3 - 3x2X3 + X2X3 + 2x2X3) 

111 



(I-X1X2) (1-X2)3 (1-Xi)5 

+(x^x- + 3x^xJ^ + 3x^xJ° + 4x1) (1-xW 

+ ( — X2xJ^ — X2X3^ — 6X2X3^ — 6X2X3° — 10X2X3° — 10X2x| — 10X2x| — IOX2X3 — 6X2X3 
-6X2X3 - X2X3 - X2X3) 

1 11 

(1 — X1X2X3) (1 — X2) (1 — x^y 

+(6x2X3 - 3xi - 3x2X3 + X3 + X2x| - 3x2X3 + 3x2x| + 2xix| - 2xix3 + 4x2X3 
-x|x|) 

111 



(1 - Xl) (1 - X2X3) (1 - X3)'' 

+ (-3X2X3 + 2X2X3 + 2X2X3 - X2X| - X2X|) 

+(-2X2X3 + X2X3 + X2X3) 



(1 — X1X2X3) (1 — X2X3) (1 — X3)^ 
1 1 



(1 — X1X2X3) (1 — X2X3)2 (1 — Xs)*^ 



2 4 1 1 1 

-X2X3- 



' (1 — X1X2X3) (1 — X2X3)^ (1 — Xs)^ 

+(x2xj° + X2X3 + 6xix| + 6x2x| + 10x2x| + 10x2X3 + 10x|x| + 10x2x| + 6x|x3 + 5x2X3 

+xlxl + X2X3) 

1 11 

(1 - X1X2X3) (1 - X2Y (1 - xlY 



+ (— X2X3 + 2xix3 — 2x\x% + 3x2xf — x| — 4x2X3 + 9xix3 — 9xix| + llx|x| — 6x3 
— 6x2x| + 15x2X3 — 15x2x1 + 15x2X3 — 9x3 — 4x2X3 + 11x2X3 — 11x2X3 + 9X2X3 

r, 2 2 I o 2 3 ,5 3 5 , n 5 4 r, 4 10 , 6 8 4 9 , o 6 7 

—7X3 — a;'2X3 + 3x2X3 — 3x2X3 + 2x2X3 — 2x3 — X2X3 + X2X3 — 3x2X3 + 3x2X3 
-3x|x| + 3xix| - X2x| + xixi) 
111 



(1-Xi) (1-X2)3 (1-Xi)5 

+ (12X2x| + X2X3 + 3x2X3''' + 10X2X3 + 6x2X3) 
+ (3X2X3 - 3X2xi - 2x1x3 + X2X3 + X2X|) 



1 11 

(1 — X1X2X3) (1 — X2) (1 — X3)''' 

111 



(1 — X1X2) (1 — X2X3) (1 — X3)'' 



, / 3 4 2 3x 1 1 1 

+ (X2X3 - X2X3)- 



(1 - X1X2) (1 - X2X3)2 (1 - X3)'^ 



+ ( — 3X2X3 — X2X3 — X2x| — 3X2X3)— — 



'(1-X1X2X2) (1-X2)2 (1-Xi)6 
+ (— X2X3^ — 6X2X3^ — 16X2X3^ — 20X2X3" — 16X2x| — 6X2x| — X2X3) 



(1 — X1X2X3) (1 — X2X3) (1 — X3y 



+ ( — 2X3^ — X2X3^ — 11X3^ — 4X2X3^ — 21X3" — 6X2X3"*^ — 17X3 — 4X2x11' — 5X3 — X2X3 + 6X2X3""^ 

, 1 ™ 3 12 r, 2 13 I n 3 15 , , 2 12 , , , 3 14 01 2 11 , 01 3 13 i ^ 2 10 
+ 17X2X3 —2X2X3 +2X2X3 — 11X2X3 +11X2X3 —21X2X3 +21X2X3 —17X2X3 

r 2 9 I 4 16 , o 4 15 , o 4 14 , 4 13\ 
— 5X2X3 + X2X3 + 3X2X3 + 3X2X3 + X2X3 ) 

1 1 1 



(1 - Xl) (1 - X2X2) (1 - xjy 

+ (10x2X3^^ + 50x2x" + 105x2X3'^ + llQxlxf + 77x2x" + 27x2X3° + 4x2x| + IQxlxf 
+40xix" + 63xixf + 49x|xJ^ + 19xix" + 3x1x3° + ^xlxf + 17x1x3"* + 21xlxf 

+ 11X2X3^ + 2X2X" + X2°X3® + 3X2°X3'* + 3X2°x" + X2°X3^) 



(1 — Xix|x|) (1 — X2X3) (1 — x|) 



2U 
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+{-2x1x1'' - llxlxl'^ - 2\x\xf - Vlx\xf - ^x\xf - x\xf - Zxlx"^ - 2,x%xf 

6 15\ 
-0:2X3 ) 

1 11 

(1 - X.xlx^^ (1 - ^2^2) (1 _ ^2)7 

+(— X2^X3* — 4x2"x3'* — 6x2X3" — 4x1x3* " X2X3'* — Gx^'xs^ — 14x2X3^ — 16x|x" 

n 7 13 r> 6 13 n 9 12 8 12 -, . 7 12 ^ 6 12 5 12 11 13 

—9X2X3 — 2X2X3 — 9X2X3 — 16X2X3 — 14X2X3 — 6X2X3 — X2X3 — X2 X3 

— 2X2°X3^ — X2X3^ — 4X2X3^^ — 6X2X3^ — 4X2X3^ — X2X3^) 
111 

(1-xi) (i-xixi)3 (i-xi)5 

, /n 3 13 2 12 7 , 1 1 3 12 ^ 2 11 10 , 01 3 11 e 2 10 o 9 

+ (2X2X3 — X2X3 —X2X3 + 11X2X3 —4X2X3 — X2X3 +21X2X3 —6X2X3 -8X2X3 

+ 17X2X3° - 4X2X| - 8X2X3 + 5X2x1 - X2X3) 
1 1 1 

(1 — X1X2) (1 — X2X3) (1 — xl)'^ 

+{xlxf + X2xJ^ + 5xixJ^ + 5x2xJ^ + lOxix^^ + 10x2xJ^ + lOx^xJ^ + 10x2xJ° + 5x^xJ' 

+5X2x1 + X2xl + X2X3) 

1 11 

(1 — X1X2X3) (1 — X2X§) (1 — xl)'' 

+ (X2X3* + 6X2X3^ + 12X2X3^ + 10X2X" + 8X2xJ°)--; ^- 5^- 7- ^- ^7- ^-^^ 

(1 — X1X2X3) (1 — X2X|) (1 — X|)' 

+(-3x|xi + 2xix| - 4x1x1 + -1-^(13^^(13^^(1^ 

+ (-xix| - Xixl) ^^——^^ (1-X2X3)2 (1 - X3)« 



+ (X2^X3^ + 4X2^X3^ + 6X2^X3® + 4X2°X3® + X2X3® + 6X2^X3® + 14x"x3® + 16X2°X3® 

, n 9 15 , „ 8 IS , „ 11 14 , , „ 10 14 , 1 , 9 14 , „ 8 14 , 7 14 , 13 16 
+9X2X3 +2X2X3 +9X2 X3 +16X2 X3 +14X2X3 +6X2X3 + X2X3 + X2 X3 

+2X2^X3" + x"x" + 4X2°X" + 6X2X3^ + 4xix3^ + X2X3^) 

1 1 1 

(1-Xlxixi) (1-Xixi)3 (1-Xi)5 

+ (X2xi - X2X3) 



(I-.T1) (1-X2)4 (1-X3)4 

+(— 3xixi + 8x2X3 + 2x2X3 — X2X3 — X2X3) 



(1 - X1X2) (1 - X2) (1 - Xs)'' 

+(2xixi — 2x2X3 — 2xixi + X2X3 + X2X3)- 



32r,2 r,33,22 ^ 1 1 1 



(I-X1X2) (1-X2)2 (1-X3)« 



4 



, / 6 5,6 6n 1 1 1 

+ (-X2X3 + X2X3) ^^-^-^ (1 - X2)4 (13^3) 

, /o 5 7 „ 5 8 , . 6 8 6 9n 1 1 1 

+ (3X2X3 - 2X2X3 + 4X2X3 - X2X3) ^^-^^ ^^-^ ^^-^ 

+ (-X2X3 + X2 + 2X2X3 - X2X3 - X2' ""^ 



(1 — X1X2) (1 — X2)3 (1 — Xs)^ 

+(3xixi - 2x1x3 - 2x2x1 + x\x\ + X2X|) 



4s 1 11 



(1 — X1X2X3) (1 - X2) (1 — Xs)'^ 



,2223 224 'i\ 1 1 1 

+ (— X2X3 + X2X3 + X2X3 — X2X3 — X2X3) 



+ (8X2X3 - X3 - X2X3 - X2X3 + X2X3 + X2X3) 



(1 - X1X2X3) (1 - XlY (1 - X3)'^ 
1 1 1 



(1 - Xl) (1 - XlXzY (1 - X3) 

^31 1 1 

+X2X3 



(1 — Xl) (1 — X2X3)3 (1 — X3)^ 

+(-6x2X3 + 8x3 + 3x2x1 - xl - X2X3 + 3xixl - 3xix| - 2x1x3 + 2x^x1 - 4x|xl 

+x|x|) 

111 

(1-Xl) (I-X2) (1-X3)7 
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+(3a;2a;3 — 2x3 — 2x2X3 + X3 + X2x\ — 2x\x% + 2x\x\ + 2x\x\ — 2x\x\ + 2>x\x% — x\x^^ 
111 



(1-Xl) (1-X2)2 (1-X3)« 

+{—X2 + xlxs — xlxl + X2X3 — 2x1x1 + 2x1x1 + 1 — 2xtxt + xtxl) 



+ (X2X3 — 2X2X3 + 2X2X3 — X2X3) 



(1-Xi) (1-X2)3 (1-X3)S 
1 1 1 



(l-xixlxl) (1-0:2)3 (l-xa)^ 
+{-2xlxt + 2x1x1 -^xtxl+xtxl)- ^ ^ ^ 



'{l-x-LX^xD (l-a;2)2 (l-a;3)6 



TT - — - — = +(—xtxl^X3 — lOxf a;2^a;| — 23xtxl°X3 — 2x\x^x% — 8xia;2^a;| — lQx\x^x% 

I — x" 

aeC3 

4 12 8 r> 3 11 8 2 10 8 < 6 11 8 orv 5 10 8 4 9 8 

—X1X2 X3 — 2x^X2 X3 — X1X2 X3 — 4X1X2 X3 — 20X1X2 X3 — 62X1X2X3 

on 9 8 . 2 9 8 ^ 6 10 8 on 5 9 8 00 4 8 8 on 3 8 8 

— 2UX1X2X3 — 4X1X2X3 — DX1X2 X3 — 2UXXX2X3 — 20X^X2X3 — 2UX1X2X3 

^288 ,6 9 8 io588 o478 10 3 78 ,278 

— DXXX2X3 — 4X1X2X3 — IUX1X2X3 — 8X1X2X3 — IUX1X2X3 — 4X1X2X3 

6 8 8 ,4 6 8 o 3 6 8 , 2 6 ,8 n .6 ,10, 6 

— X1X2X3 — ZX1X2X3 — X1X2X3 — ZX1X2X3 — X1X2X3 — 0X1X2 X3 

-18x1X2X3 - 30xfx2x| - 6xiX2°X3 - 18x^x|x| - 18x?xix3 - 3x^X2°x| 
— 6x?xix| — 3x?x|x| — &x\x%x% — 18xf a;|a;| — 18a;iX2x| — 18x?a;2a;| 

—6x1X2X3 — 3x1X2X3 — 6x^x2X3 — 3x1X2X3 — 6x1X2X3 — 3x^x2X3 

, o 6 11 7 , o/i 5 10 7 I .0 4 9 7 I ^ 5 11 7 , 01 4 10 7 I o,f 3 9 7 
+3X1X2 X3 + 24xiX2 X3 + 48X1X2X3 + 6X1X2 X3 + 21X1X2 X3 + 24x1X2X3 

+3XiX2"'^X3 + 6X1X2*^X3 + 3X1X2X3 + 9X1X2*^X3 + 36x^X2X3 + 48xix|x3 

+36x?a;2a;3 + 9x1X2X3 + 9x1X2X3 + 24xf X2X3 + 21x^X2X3 + 24x^x^x1 

+9X^X2X3 + 3X^X2X3 + 6X1X2X3 + 3X1X2X3 + 6x?X2X3 + 3X1X2X3) 

1 11 

(1 — x\x\x3) (1 — X^Xs) (1 — X3)'^ 

+(3x2''x| + 18x2^x| + 45x2^x| + 60x2^x| + 45x2''x| + 18x2x| + 3x2x| + 10x2^X3 
+40x"xl + 60x2''x;^ + 40x1x3 + lOxIxg - 10x2^x| - 50x2^x1 - 100x"xi - 100X2°x| 
-50x1x1 - 10x^x1) 

1 11 

(1 — X1X|X3) (1 — X^Xs) (1 — X3)'' 



,/ 98,, 88, c 78,. 68, 58 10 8 £-98 1088 

+ (XlX2X3 + 4X1X2X3 + 6X1X2X3 + 4X1X2X3 + X1X2X3 — X2 X3 — 5X2X3 — 10X2X3 

— 10x2x| — 5xlx| — X2xi) 

1 11 

(1 - X1X2X3) (1 - X|X3) (1 - Xs)'^ 

+(2x1X2X3 — 2xiX2x| + 11x1X2X3 — 11x1X2X3 + 24x1X2X3 — 24x\x%x\ 
+26xix|x| — 2%x\x%x% + 14x?xlx| — \4x\x\x% + 3x?x|xi — 3xiX2xi 

,396,<386,£.375,.365, 355 11 6 £.10 6 
+X1X2X3 + 4X1X2X3 + 6X1X2X3 + 4X1X2X3 + X1X2X3 — X2 X3 — 6X2 X3 

— 15xix| — 20x1x3 — 15x2xf — 6xix| — xix| + x|xi + 5x2x1 + 10x|x| 
+10X2X3 + 5X2X3 + X2X3) 
1 1 1 

(1 — X1X2) (1 — X|X3) (1 — X3)'' 

+ (-X?X2xi - 3x?xixi + 3X1x1x1 - x\x\ + x|x| + 5xixi) 7- ^ r — ^- r — ^- 

(1 — X1X2) (1 — X2X3) (1 — X3 

I / 4 7 5 I o 3 6 5 I 2 5 5 . q 4 6 5 . ^ ?> 5 5 , q 2 4 5 
-\-\X\X2X3 -j- ZX1X2X3 -t" X1X2X3 + 0X1X2X3 + 0X1X2X3 + 0X1X2X3 

+3x^X2X3 + 6x?X2X3 + 3x\x\x\ + x\x\x\ + 2x\x\x% + x\x%x\ 

—2x2X3 — 10x2X3 — 20x2X3 — 20x2X3 — 10x2X3 — 2x2X3 + X1X2X3 + 3xiX2xi 

+3X1X2X3 + XlX2x| — X1X2X3 + X1X2X3 — 5X1X2X3 + 5X1X2X3 — llx\x\x\ 
+9xixix3 — UXixixs + 7X1X2X3 — 4Xixix3 + 2X1X2X3 — X1X2X3 — 7xfx2X3 

— Ilx?x|x3 — 5x?xix3 — 2xtx%x\ — 4xixix3 — 2x\xtx\ — x%xt 

—4x2X3 — 6xix3 — 4x1x3 — X2X3 + 33xlx| + 22x1x1 + 3x|xi — 3xlx| — xlx| 

+5x1x1 + 21x1x1) 

111 

(1-Xl)2 (1-Xix3) (1-X3)« 
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+ {xlxl'^xl + Qx\xlxl + Qxixlxl + xlxl'^xl + 2x1x1x1 + x\xlxl 
+Axlxlxl + lAxtx^xl + 16x1x1x1 + Axfxlxl + 6x1x^x1 + IQxjxlxl 
+14xlx2X^ + 6xiX2xl + Ax\x\x\ + 9x\x%x\ + 6x\x\x\ + Ax'ix\x\ 

,667,r>467, 347, 247 ol2 8 10 H 8 ^rlOS 
+X1X2X3 + 2X1X2X3 + X1X2X3 + X1X2X3 — 6X2 X3 — 19,X2 X^ — 45X2 X2, 

—QOx%x%, — 45a;f a;| — 182:22^3 — Sxf^^l — xix%x\ — Ax\x\x\ — 6x1X2X3 — 4x1X2X3 
—xix^xl — 2x1X2X3 + 2x1X2X3 — lla;?a;2a;3 + lla;ia;2a;3 — 21x\x2x\ 
+24a;ia;2a;3 — 2Qx\x2X3 + 26x\xtx'3 — lla;ia;2a;3 + 14:X\X2X3 — 'ix\x\x% 
+ixix'ix% — x\x%x% + 2x1x1x1 + 9x1x2x1 + 8x1x1x1 + 2xlxtx3 
+3a;i.T2.T3 + 12xixlxl + 3x^X2X3 + 6xixlxl + 15a;ia;2a;3 + 3x^X2X3 

—3xixlx3 — 9x1X^X3 — 27x\x\x% — 27x'ix\x% — 9x\x%x% — 9x'{x\x% 
—21x\x%x% — \^x\x%x% — 9x\x\x% — 3x\x%x% — 6x1X2X3 — 3x1x2X3 
—3x\x\xf^ — .T2.T3 — 5.T2X3 — IOX2X3 — lQx2x'3 — 5a;2a;3 — x%x\ + 101a;23;3 
+105x^3:1 + 6Qxlxl + 2Qx%xl + hxlxl + xWi + 10x"xJ + 50a;2°a;I - 9x^^x1 
-34x1x1 - 45iM - 20.4;4 + ^xlxl + 6x1x1 + xtxl) 

1 1 i 

(l-a;i) {l-xlx3) {l-X3y 

+{x\x\xl + 2x\x%x% + x\x%x% + 2x1x1x1 + Axfx^xl + 2x1x^x1 

,478,r,368, 26 8n 
+X1X2X3 + 2X1X2X3 + X1X2X3) 

1 1 1 

(1 — x\x\x3) (1 — X2X3) (1 — X3y 

+{-'^xlx\\l - 12a;?x^xi - 2Qx\xlxl - 4x\xfxl - \2x\xix% - 12x?x^a;i 
—2xixi''x3 — 4x^X2X3 — 2x1X2X3 — 4x^X2X3 — 12a;fa;2a;3 — 12x\x2X% 
—12x1x1x3 — 4xixlx3 — 2x1X2X3 — 4x\x\x% — 2xtx%x% — 4x\x%x% 
—2x'ix%x% + xtxl^xl + 8xlxl°xl + 16xixlxl + 2x1x1^x1 + 7xtxl°xl 
+8xlxtxl + xixi^xl + 2xlxl°xl + xlxtxl + 3xtxl°xl + 12x1x1x1 
+16a;ia;2a;3 + 12a;f a;2a;3 + 3a;ia;2a;3 + 3xiX2xl + 8x^X2x1 + 7x\x''2x\ 
+8x\xlx''3 + 3x1x1x1 + x\x%xl + 2x1x1x1 + x\x%xl + 2xixlxl 
+xixlxl) 

1 11 

(1 - Xlxlx3) (1 - X^X3)^ (1 - 0:3)'^ 

+{—Xixl''x3 — 6a;f a;|a;| — 10a;ia;|a;| — 2xixl°X3 — 6a;ta;|a;| — 6a;fa;|a;| 

4 10 6 n396 286 r,696 c586 r;476 
—X1X2 X3 — 2X1X20:3 — X1X2X3 — 2X1X2X3 — 6X1X2X3 — 6X1X2X3 

c376 o276 686 ^576 466 r>366 
—6X1X2X3 — 2X1X2X3 — X1X2X3 — 2X1X2X3 — X1X2X3 — 2X1X2X3 

-xJxaXg) 

1 11 

(1 - Xfxixs) (1 - Xix3)3 (1 - X3)5 
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,/ 684, 4 5 4,o464 o584 o474,r,354 
+ {—XiX2Xj + X1X.2X1:. + ,iXiX2X^ — 2XiX2X-^ — ,iXiX2X-^ + 1X1X2X3 

-xjxlxs - 2x'ixlx\ - x\x%x\ + xlxlxt + 2x\x%xl + xix\xl) 

1 11 

{l-xlxlx,,) {I-X2Y {l-xsY 

+ {x\x:2X'i, + AxiX%x"i, + 4X1X2X3 + 2xiX'2X'^ + + 4X1X2X3 

+xixlxl + 2x1x1x3 + xixlx3 + Xixlx's — xtx2X3 — 3xia:2a;3 

r>366, 266 676 r>665 255 684 

— 2a;ia;2a;3 + X1X2X3 — X1X2X3 — 2x^X2X3 — X1X2X3 — X-1X2X3 

r,5 74 464 o684 <474 r,364 484 

— 2X1X2X3 — X1X2X3 — 2x^X2X3 — 4X1X2X3 — 2X1X2X3 — X1X2X3 

-2xixlxt - xixlxt) 

1 11 
(1 - xfxixs) (l-a;2)3 (l-xsY 

, I 6 10 6 /;596 in486 r,5 10 6 ^49 6 ^,3 86 
+ (— X1X2 X3 — bXxX2X3 — \yiXxX2X3 — 2xxX2 X3 — bXxX2X3 — bXiX2X3 

4 10 6 r, 3 9 6 2 8 6 o r, .6 .4 J,fi 

— X\X2 X3 — AX\X2X3 — X ^X2X3 — AX-^X2X3 — V)X-^X2X3 — 0X^X2X3 

£.376 o27 6 6 86, 456,o 4 6 6|o356 

— oa;iX2X3 — 2x1X2X3 — X1X2X3 + X1X2X3 + 6x1X2X3 + 2x1X2X3 

266, 676, 0666, 266 684 r>574 
—X1X2X3 + X1X2X3 + 2X1X2X3 + X1X2X3 — X1X2X3 — 2X1X2X3 

-Xi.x-2.-4 - 2x\x%x\ - 4x^x1x3 - 2x\xlxl - x\xlxl - 2x\xlx\ 

— X1X2X3 + 2x^X2X3 + %x"lx%X3 + lQx\x\x"3 + 4x"{x2X"3 + lQx\x%x\ 

+d,xixlxl + 2xtxlx'3 + 4xixlxl + 2x1x1x1 + 2x\x\xl + 4x\x\xl 
+2xtxtxl + 4x1x1x1 + 2x'ix%x\) 

1 1 1 

{1-Xlxlx3) (1-0:2)2 (1-0:3)6 



,/6 117,o5 10 7,ic-197,r,5 117,~4ia7,„397 

+ (a:ia;2 X3 + 8.X-1.X-2 X3 + lDa;iX2X3 + 2x1X2 X3 + 7x1X2 X3 + 8x1X2X3 

, 4 11 7 I n 3 10 7 , 2 9 7 , o 6 10 7 , , „ 5 9 7 , , ^ 4 8 7 
+X1X2 X3 + 2X1X2 X3 + X1X2X3 + 3XiX2 X3 + 12x1X2X3 + 16X1X2X3 

+12a:ia:ia:3 + 3a:?a;|a;3 + 3a;fa;2a:3 + Sxfalxs + Gxfxlxl + 8a;?a;2a;3 

,o277, 687 4 5 7 r,467 o357, 267 
+3X1X2X3 + X1X2X3 - X1X2X3 - 3X1X2X3 - 2X1X2X3 + X1X2X3 

X1X2X-3 ^X 1 .-J/-2 .'^-,3 0:1X2X3 tJX-iX2X3 ^^X 1X2X3 'X' 10X1X2X3 

+ 6X1X2X3 + 15X1X2X3 + 12X1X2X3 + 3X1X2X3 + 6X1X2X3 + 3X1X2X3 

+3X^X2X3 + 6X^X2X3 + 3X1X2X3 + 6X^X2X3 + 3x^X2X3 — 3XiX2°x| 

-18x?x|x| - 30X1X2X3 - 6x?X2"x| - ISxiX^X^ - 18x?x|x| - 3XiX2°X3 

— 6x?x|x| — 3xfx|x3 — 6xiX2x| — ISx^x^xl — 18X1X2X3 — 18x?X2x| 
— 6xfx2xf — 3xf x|x3 — 6xf X2X3 — 3xtx2xf — 6x?X2xf — 3x^X2X3) 

1 1 1 



(1 - X^X^X3) (1 - X2) (1 - X3) 

+{xlxl - xtxl) ^ 



7 



(1 — X1X2X3) (1 — X2)® (1 — Xs)^ 

1 1 1 



+ ( a^2+0:2)-(^_^^^^)(^_^^^6(^_^^)2 
+ ( — 2X2X3 — 2X2X3 — X1X2 + X1X2 + X?X2 — Xl — 2X2 + 1 + 2X2X3 + X2X3 + 2X2X3 + X2X3 — X2X3) 

1 1 1 

(1-Xl)2 (1-X2)4 (1-X3)3 
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+ {3xlx'l + Sx2xl + 2x\x\ — 2x\x\ — 2x\x% + 2xiX2 + 2x1 — X2 — 2xtx\ + x\x\ — 2x%x% 

-ix\xi - X^Xi - x\x3 + X2X3) 
111 



3 



(l-Xl) {l-X-zf {I-X3) 

+{-hx%xl - 3a;2a;|) ^ 



(1 — XlX%X3) (1 — X2X3) (1 — xzY 

1 11 



^'""^ (1 - xixlxs) (1 - 0:2x3)2 (1 - 0:3)6 
+(4a;2^X3 + 16a;2^a;3 + 2Ax^x\ + 16a;2a;3 + 4a;2a;3 — 2a;2^a;| — 10a;2^a;3 — 20a;2'^a;3 
-20a;2°a;i - lQxlx% - 2xlx%) 
1 1 1 

(1 - xixlx-i) (1 - xlxs)'^ (1 - X3)*' 

, / 12 7 I . 11 7 , c 10 7 I < 9 7 I 8 7\ 

+(o,-2 X3 + 4a;2 X3 + 6x2 xz + 4a;2a;3 + a;2a;3)- 



(1 - xix^xs) (1 - x^xs)^ (1 - 0:3)^ 
HSxtxl - 344) (i_ ^'^^2,3) (i_',^)e (i_',3)2 

+(-30.14 - 344 + 44 + 30.2^4 + 244) (i_^'^2^3) (TZ^ (TZ^ 

+(3a;2°a;3 + ea^Xs + 3x^x1 - 2x^x1 - 3x1x1 + x^xl - Ax%xt - 4x24) 
1 1 1 

(l-a;ixia;3) (1-0:2)'' (1-0:3)* 

+(— 3o:2^o:| — 9o:2°a;3 — Oxfxl — 3a;|o:| + 3o:2o:| + 3o:2a;| — 3a;|a;3 — 30:20:3 
+6o:^°o:i + 12o:lo:i + 6o:^o:i) 
1 11 



(1 - 0:10:^0:3) (1 - 0:2)3 (1 - 0:3)5 



+ (3X2^0:3 + 12x"x3 + 18X2''X3 + 12X2X3 + 8X2X3 - 40.-2X3 - 8X2X3 + 6X2°x| 

+12x1x1 + 6xixi - 8x"xl - 24x2°xi - 24x1x1 - 8x1x1) 
1 1 1 

(1-Xixix3) (1-X2)2 (1-X3)6 

+(—3x2^x1 — 15x2^x1 — 30x"x| — 30x2°xl — 15xlx| — 3x|x| + 6x2X3 + 3x2x1 
-10x"xl - 30x2°xl - 30x1x1 - 10x1x1 + lOxa^Xg + 40x"xl + 60x2°xl + 40x1x3 
+10x1x1) 

1 11 
(1 — X1X2X3) (1 — X2) (1 — X3)'' 

, r A 4,r,2 4 435 o325 2 5 5,324 
+ (— 4X1X3 + 2.X-1.X-2X3 — X1X2X3 — 2X1X2X3 — X1X2X3 — X1X2X3 +X1X2X3 

+4x3 + X2xl + 2x1x1 + 3x1x1 - 4x2X3) 

111 

(1-Xl)2 (I-.T2X3) {l-Xzf 

+ (8X1X2X3 + 3X1.T2X3 + .T1X2X3 — X^X2X3 — 2XiX2x| — x'lXjjxJ 

457 r,3 47 2 3 7, 37 r 5 37 a 2 6 

-X1X2X3 - 2X1X2X3 - X1X2X3 + X1X2X3 - 5X2X3 - X2X3 - 4X2X3 

+8X2X3 + 5X2X3 - Xjxl - 5x2x1) 

111 

(1 - Xl) (1 - X2X3) (1 - X3)'^ 

+ (-3xixixl - 3xixlx3 - Xix|xl + 6x1x1 + 4x1x3 + x\x\)— ^ r--- ^ r--- ^ ^ 

(1 — X1X2X3) (1 — X2X3) (1 — X3)^ 

+(-2x1x1x1 - xixlxl + 3x1x1 + xlxl) (^3^^ (l-xU)^ (T^ 
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(1 — X1X2X3) (1 — X2X3)^ (1 — Xa)^ 

1 



+{xiX2Xs - X2X3 - X2X3 - X1X2X3 + X2X3 + X2X3) 



(1 - X1X2X3) (1 - X2y (1 - 2:3) 



4 



,/• B5 55,75,r>D5,55, 34, 45 2a 

+ {-XlX2Xs - X1X2XS + X2X3 + 2X2X3 + X2X3 + X1X2X3 + X1X2X3 - X2X3 

-3 
1 



-2xlxt - xtxl) 



(1 - X1X2X3) (1 - X2)^ (1 - Xs)^ 

1/ 76, o 66, 56 86 o76 066 56 24 

+(a;ia;2a;3 + 2a;ia;2a;3 + a;ia;2a;3 - a;2a;3 - 3a;2a;3 - 3a;2a;3 - a;2a;3 - X1X2X3 
—2xixlxl — xixtxl + 3xlxi + 3x1x1 + xtxl) 
1 11 



(1 — X1X2X3) (1 — X2)^ (1 — X3)^ 

, I 87 o 77 o fi7 5 7,97,<87,c77,<67 

+ y-XxX2X3 - ZX\X2X3 - Zx-iX^X-,, - X1X2X3 + 0:2X3 + 40:2X3 + 6X2X3 + 4^2X3 

+X2X3 + 3X1X2X3 + 3X1X2X3 + X1X2X3 - 6X2X| - 4X2X3 - X2X3) 

1 11 

(1 — X1X2X3) (1 — X2) (1 — X3)'' 

+ (X1X2X| + 2X1X2X3 - X2X3 - 3X2X| + X2X3 - xlxl)j:^ 

+ (X1X2X3 — X2X3) 



(1-Xl) (1-X2X3)2 (1-X3)« 

1 1 



(1 — Xl) (1 — X2X3)^ (1 — X3' 



5 



+(-XiX2X3 + X1X2X3 - 4xixlxj + 4xix|x3 - 6x^X2X3 + 6xix|x3 - 4xix|x3 

+4xix|x3 - xlx\x\ + X1X2X3) 

111 

(1 - XlX2)2 (1 - X|X3) (1 - X3)'^ 



+ (XiX2X3 — X1X2X3 + 3X1X2X3 — 3X1X2X3 + 3X1X2X3 — 3xiX2x| + x\x\x3 

3 3 2 3 3 I 2 3\ 

-X1X2X3 - X1X2X3 + X1X2X3) 

1 1 1 

(1-X1X2)2 (I-X2) (1-X3)6 

+(-2xlxt + 2xixl + xtx3 + X2X3 + 2x?xi - 2xiX2 - 2x1 +x%- xtx3)j- — rj- — - — rrj- — - — r 

(1 - X1X2) (1 - X2)^ (1 - Xs) 

+ (2x^X2X3 — 2X1X2X3 + 3x?x|x3 — 3X1X2X3 + x\x\x3 — x\x\x3 — 3x?X2X3 + 3X\X%X3 

72r,62 52,4 ,o3 2 ,4 2-, 

— X2X3 — 2X2X3 — X2X3 + X2X3 + 3X2X3 — X2X3 + X2X3) 

1 1 1 

(I-X1X2) (1-X2)4 {I-X3Y 

,/o262,r, 52 r252,r 42 o242,o 3 2 362 
+ ( — 2X1X2X3 + 2X1X2X3 — 5X1X2X3 + 0X1X2X3 — 3X1X2X3 + 3X1X2X3 — X1X2X3 

352, 342,o232 o 22, 83,o73|o63 
— X1X2X3 + X1X2X3 + 3X1X2X3 — 3X1X2X3 + X2X3 + 3X2X3 + 3X2X3 

+X2X3 - X2X3 - 2X2X3 - 4x1x3 + xlx3 - X2X3) 

111 

(1 — X1X2) (1 — X2)^ (1 — Xs)^ 

+{2x\x\x% - 2xix%x% + 7x1X2X3 - 7xiX2x| + 8xixlx| - 8x1X2X3 + 3xiX2x| 

— 3X1X2X3 ~l~ X1X2X3 ~l~ 2X1X2X3 ~t~ X1X2X3 — X1X2X3 — 3X1X2X3 

+3xixlx3 — X2X3 — 4x1x3 — 6x2X3 — 4x1x3 — xlx3 + xlx3 + 3x1x3 + 3xlx| 
+5x1x1 - xixl + xlxs) 

111 

(I-X1X2) (1-X2)2 (1-X3)« 
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+ {—2xiX2xf + 2xixlx3 — 9xixlx3 + 9x1X2X3 — lZx\x%x'k, + l^X\X2x\ — 11x1X2X3 
+lla;ia;|a;3 — 3xiX2xt + ^xix^x^ — XiX2xt — Sxfxlxt — 3x1X2X3 
-xlxlxt + xlxtxt + 3xixlxt - 3xixlxt + xl°xl + 5x1x1 + lOx^xl 
+10x1x3 + 5x^X3 + X2X3 — x\x\ — Ax\x\ — &x\x3 — Ax\x3 — 6x%x\ + x'ixt 
-xtxl) 

1 1 1 

(1 - X1X2) (1 - X2) (1 - XsY 

, I 262, 52 „262,r, 42 242, 32,232 

-\-\~XxX2X3 + X\X2X3 - 2X1X2X3 + 2X1X2X3 - X1X2X3 + X1X2X3 + X1X2X3 

2 2n 
-X1X2X3) 

1 1 1 



(1-X1X2)2 (1-X2)2 (l-Xa)^ 

+{x'ixlx3 — Xixtx3 + xixlx3 — X1X2X3 — xfalsa + Xix'ix3) 



{I-X1X2Y (1-X2)3 {I-X3Y 

24. ;^.2;s 2\ 1 1 1 

~r X'iX'2 ~r X1X2 X1X2) 



'(l-xiXiJ^ (1- 2:2)4 (1-2:3)3 
+{-xtxtxt - 2xlxlxt - x'ixtxt - 2x\xlxl - 'lx\x\xl - 2xixlxt 
-xtxtxt - 2xlxlxt - x\xixl + a;2a;3 + Ax^xl + &X2x\ + Ax\x% 
+x'lx\) 



(l-XlY {1-XIX3Y {I-X3Y 

,/ 442 r»332 222,o63,,f53,o43 22 

+{-XiX2X3 - 2x1X2X3 - X1X2X3 + 2a;2a;3 + 40:2X3 + 20:2X3 - 0:10:20:3 

+x\x%X3 — X1X2X3 + 2x\x\x3 — 2X1X20:3 — 2X1X2X3 + 2X1X3 + X1X2X3 + X1X2X3 + 2X1X2X3 
+X1X2X3 + X1X2X3 — X?X2X3 + X2X3 + 3X2X3 — 2X3 — 2X2X3 — 4X2X3 — X2X3 — 3X2X3) 



(1-Xl)2 (1-X2)3 (1-X3)4 



,/-453,r,343, 233, 443,r,3 3 3, 22 3 
+ (XiX2X3 + 2X1X2X3 + X1X2X3 + X1X2X3 + 2X1X2X3 + X1X2X3 

n74 e64 c 5 4 o**, 33, 23 „ 2 242 

— 2X2X3 — 6X2X3 — 6X2X3 — 2X2X3 + X1X2X3 + X1X2X3 — 3X1X3 — X1X2X3 

I 32 o232io 22 o222if^ 2ii^2 2 342 

+X1X2X3 — 3X1X2X3 + 3x1x2X3 — 3x1X2X3 + 2x1x2x3 + 2x1x2x3 — X1X2X3 

o 3 3 2 4 3 3 n 3 2 3 2 3 3,3 2 2 4 4 2 

— 3X1X2X3 — X1X2X3 — 2X1X2X3 — X1X2X3 — X1X2X3 + X1X2X3 — X1X2X3 

-X2X3 - 2X2X3 - 4X2X3 + 3X3 + 5X2X3 - X2xi + X2X3 + X2X3 + 4X2X3 + 9X2X3 
+2X2X3) 

111 
(1-Xl)2 (1-X2)2 (1-X3)5 

+(-XiX2xt - 2x^X2X3 - X1X2X3 - 2xix|x3 - 4x?X2X3 - 2xlxlx\ 
— X1X2X3 — 2xfxix3 — X1X2X3 + 2x|x| + 8x2xf + 12x%x\ + 8x|xf 

,o4 5 44 „ 34 24,. 3,2 5 3 43, .243 

+2X2X3 — X1X2X3 — 2X1X2X3 — X1X2X3 + 4X1X3 + X1X2X3 — X1X2X3 + 4X1X2X3 

. 3 3,^ 2 3 3 r 2 3 , < 2 2 3 r, 3 „ 2 3,3 5 3 

—4X1X2X3 + 7X1X2X3 — 0X1X2X3 + 4X1X2X3 — 2X1X2X3 — 2X1X2X3 + X1X2X3 

+6X^X2X3 + 5X^X2X3 + X1X2X3 + 2x\x\x\ + X1X2X3 + X1X2X3 — X?X2X3 

+2X1X2X3 + 2X1X2X3 + X2X3 + 3X2X3 + 3X2X3 + 5X2X3 — 4X3 — 17X2X3 — 5X2X3 
+2X2X3 — 2X2X3 — X2X3 — 5x2X3 — 16x2X3 — 2x2x|) 

111 

(1-Xl)2 (I-X2) (1-X3)« 

+ (X?X^X^ - ^1^^^3) (l_,\,^)2 (l_Lx3)(l-X3)a 
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+ {2x1x^x1 + Sxjxlxl + lOx'txlxl + 2xtx2xl + Ax^xlxl + 2x\x%x% 
+4x1x1x1 + 10x^X2X3 + Sxfxlxl + Axfxlxl + 2x\x%xl + Ax\x\xl 
+2x\x\x% + 2x\x\x\ — x\x\x% — hx\x%x% — Yxlx^x^ — xtx2X3 

r) 3 8 ,6 _2_.7_6 n, 5 ,8, 6 n O, .6 •3,.2 .6 6 

— ^X\X'2X't^ — X-^X^X't^ — OX^X^X':^ — uX^X^X'^ — \JX ^X2X'^ — OX^X^X':^ 

-3x1x1x1 - 7x\xlxl - bxixlxl - "ixlxlxl - xlx^xl - 2x\xlxl 
—x\xtx% — xixtx% + 2x2^X3 + 10x2°X3 + 20x|x3 + 20x2X3 + 10x2X3 + 2x2X3 
-4x2°xl - 16x1x1 - 24x1x1 - 16x^x1 - 4x|xi) 
1 1 1 

(1-Xi) (1-X2X3)2 (1-X3)6 

+(xixixi + 4xfx2x| + 5xfx2xi + Xixixl + 2x1x1x1 + xf x|x| 
+2xf X2X3 + 5x1X2X3 + Ax\x\x\ + 2x1X2X3 + xf X2X3 + 2x\x\x% 

,345, 245 10 6 ^96 e86 a 7 6 66n 
+X1X2X3 + X1X2X3 — X2 X3 — 4X2X3 — 6X2X3 — 4X2X3 — X2X3) 

1 1 1 

(1-Xl) (1-Xix3)3 (1-X3)5 

-\-{x\x%x% + x\x\x% - x?X2xi + xtx%x% + 2x\x%x% + xix\x% 

— 3x2x| — 6X2X3 — 3x^X3 — Xix|x3 — 2Xix|x3 + 2X1X2X3 — 8x^X2X3 + 8X1X2X3 

,o22 o 34 553r,443 333 233 

+8X1X2X3 — 8X1X2X3 — X1X2X3 — X1X2X3 — 2X1X2X3 — X1X2X3 — X1X2X3 

+X1X2X3 + X?X2X3 - X2X3 - 8X2X3 + X2X3 + X2X3 + 8X2X3 - X2X3 - X2X3 + 4X2X3 

+4X2X3 + 8X2X3 + 2X2X3 + X2X3 - X2X3) 

111 

(1-Xl) (1-X2)4 (1-X3)4 

,/ 5 74 „464 o354 474 r,364 254 
+ ( — X1X2X3 — 8X1X2X3 — 8X1X2X3 — X1X2X3 — 2X1X2X3 — X1X2X3 

5 6 4 4 5 4 I 3 4 4 2 4 4 , o 9 5 , n 8 5 , n 7 5 
— X1X2X3 — X1X2X3 + X1X2X3 — X1X2X3 + 8X2X3 + 9X2X3 + 9X2X3 

+3xix| + xixixs + xix|x3 — xixixi + 2x?xix3 — 2xix|x3 + 5xix|x3 

K 32io232 o 22 o222jo 2i 352i 342 

— 0X1X2X3 + 8XiX2.X-3 — 8X1X2X3 — 8X1X2X3 + 8X1X2X3 + X1X2X3 + X1X2X3 

, 5 5, 4 , „_4_4, 4 , _3_3_4 . _2_3_4 , _3 ,4 _3_3_2 . ,5 6 3 

+2x^X2X3 + X?X2X3 + xtxlxl + 2x'iX2X3 + X1X2X3 + X2X3 + 2X2X3 

+4X2X3 + 2X2X3 — 7X2X3 — 2X2X3 — 4X2X3 + X2X3 — 6X2X3 — 12X2X3 + 2X2X3 

-3x|x3 - 8x2X3 - X2X3 + x^xi) 

111 

(1-Xl) (1-X2)3 (1-X3)5 

+(xf x|x| + Axtxlxl + 5x?x|x| + x\x%xl + 2x\x\xl + xixtxl 

+2x1X2X3 + 5x1X2X3 + 4x1x^x1 + 2x^X2X3 + x\x%x\ + x\x\x\ 

3 4 5 I 2 4 5 o 10 (' 1 o 9 6 1086 1 r, 7 6 066 
— X1X2X3 + X1X2X3 — 8X2 X3 — 12X2X3 — I0X2X3 — 12X2X3 — 8X2X3 

-X1X2X3 - 2X1X2X3 - X1X2X3 + 2xiX2x| - 2X1X2X3 + 2X1X2X3 - 7X1X2X3 

+7X1X2X3 — 7X^X2X3 + 8xiX2x| — 8X1X2X3 + 8X1X2X3 + 3x?X2xf — 2X1X2X3 

— X1X2X3 + X1X2X3 + 2X1X2X3 — X1X2X3 — 2X1X2X3 — X1X2X3 

-x\x\x\ + xixixl + x?xix| + xix|x| - 2X1X2X3 - 6X^X^X3 

— 6x?xix3 — 2x^X2X3 — 4x?X2X3 — 2x\x\x\ — 2x\x%x\ — 4x\x\x\ 

— 2x?X2X3 — 2xf x|x3 — xix3 — 8x2x1 — 3x2x1 — 5x2X3 — 2x2x1 + 25x1x3 

+11x2x1 + 8x2X3 + 5x2x1 — X2xl — 3x2x1 + Sx^xl + 24x2x1 — 5x2x1 — 8x2x1 

+3x2x1 + 4x2x1 + Xjxl - Xaxl) 

111 

(1-Xl) (1-X2)2 (1-X3)« 
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+ (-a;iX2X3 " 5x1x1x1 - lx\xlxl - x\xlxl - 2x\xlxl - xlxlxl 
—ix\x%x% — 9x^x1x3 — 9xfx2X3 — 3x^X2X3 — Sxixlx^ — 7x\x%x% 

r-366 o266 566 456, 346 246 
— 0X1X2X2, — 6X1X2X2 — X1X2X3 — X1X2X3 + X1X2X3 — X1X2X3 

+3x1^x1 + I5xl'^xl + 30a;2a;3 + 30x23^3 + 15x1x1 + 3x2X-3 + xixlx^ + 3x1X2X3 
+3xiX2X^ + xiX2X% — 3xix\x% + 2x\x\x\ — 2x\X2xi, + 9x\x%x'k, — 9x1X2X3 
+\2x\x2xi, — 15x\X2x\ + 8x1X2X3 — lla;ia;2a;3 + 3x^X2X3 — 3x1X2X3 
—3x\x2X3 — xix2X3 + x^xlxi — 3xtx2xi — 6xlx2xt — 2xixtx3 

,656,r,446, 336, 236 36 o574 

+X1X2X3 + 2X1X2X3 + X1X2X3 + X1X2X3 — X1X2X3 — 3X1X2X3 

—9x\x%X3 — 3x1X^X3 — 3x^X2X3 — 6x^X2X3 + 3a;ia;2a;| + \2x\x\x% 

+ 15x\x%x"3 + 3x\x%X3 + &x\x\x\ + 3a;'xa;'2X3 + &x"lx\x3 + 15x\x^x\ 

+\2x'ix'2X3 + Qx\x2X3 + 3a;5'a;2a;3 + 6xtx2X3 + 3a;?a;2a;3 + 3a;ia;2a:3 

+a;2a;3 + Ax\x3 + &x\x\ + 4x^X3 + 6x2X3 + 5x2X3 — 61a;2a;| — 4Ax\x% — 16x%x% 

< 6 6 c 4 6 I 3 6 . , 2 ,S i n 1" 6 /in 9 6 , ri 9 S , r>c: 8 5 

— 4a;2a;3 — 0X2X3 + X2X3 + 4./:2.x-3 — lO.r^ X3 — 4^X2X3 + 9a;2a;3 + 25a;2a;3 
+20x1x1 - 3x2x1 - 5x2x1 - xtxl + xixl) 
111 



(1-X2) {i-x3y 

+ {-X-ixt +xt+ xixl - X2xt) 



{1-XlY {1-X2X3? {I-X3Y 



+(-3a;2a;3 + 3a;2a;3 + 2^2 - 2^2)- 



(l-xi) {I-X2Y (I-X3) 



2 



+{2xtx3 — 2x1x3 — xl+ xi) 



+a;22;3 



{1-XlY {1-X2f [l-Xif 
1 1 



(1 - X1X2) (1 - X2X3)'^ (1 - X3Y 
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